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Abstract 

We consider the effect of Gaussian white noise on fast-slow dynamical systems 
with one fast and two slow variables, containing a folded- node singularity. In the 
absence of noise, these systems are known to display mixed-mode oscillations, con- 
sisting of alternating large- and small-amplitude oscillations. We quantify the effect 
of noise and obtain critical noise intensities above which the small-amplitude oscil- 
lations become hidden by fluctuations. Furthermore we prove that the noise can 
cause sample paths to jump away from so-called canard solutions with high prob- 
ability before deterministic orbits do. This early-jump mechanism can drastically 
influence the local and global dynamics of the system by changing the mixed-mode 
patterns. 
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Keywords and phrases. Singular perturbation, fast-slow system, invariant manifold, dy- 
namic bifurcation, folded node, canard, mixed-mode oscillation, random dynamical sys- 
tem, first-exit time, concentration of sample paths. 

1 Introduction 

Our main focus of study are stochastic dynamical systems with multiple time scales. 
In particular, we are going to study a special bifurcation ("a folded node") in a three- 
dimensional fast-slow stochastic differential equation (SDE) with one fast and two slow 
variables. The detailed technical discussion including all relevant definitions and precise 
statements and proofs of our results starts in Section [2] In this section we want to outline 
our motivation and state our main results in a non-technical way. There are two main 
motivations for our work: 

(Ml) We want to develop an analogue to the intricate deterministic bifurcation theory for 
random dynamical systems by linking stochastic sample-path techniques and the 
well understood deterministic theory. 
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(M2) A detailed analysis of noise effects in multi-scale stochastic systems is often crucial 
in applications; in particular, many biological systems have widely separated time 
scales and are influenced by various random effects. 

We are going to describe our two main motivations in more detail, starting with 
(M2). Complex oscillatory patterns have been discovered in many different applications. 
Chemical systems [HSl ED E] and neuronal dynamics (HU [Ml EB] provide ample examples. 
Recent work has shown |32l [63] that fast-slow systems can be used to model a wide variety 
of oscillatory patterns. A classification of local and global fast-slow "mechanisms" can be 
used to analyze each pattern. Mixed-mode oscillations (MMO) alternate between small- 
amplitude oscillations (SAOs) and large-amplitude oscillations (LAOs). Figure [l] shows 
a typical MMO time series where the deterministic time series has been generated by an 
MMO model proposed in [29]. The time series shows alternations between L = 1 LAOs 
and s = 7 SAOs which is denoted as the MMO pattern = . 




Figure 1: Comparison of a deterministic (black) and stochastic (red) time trajectory. 
The left-half of the figure shows the time series exhibiting a 1^ MMO; see also the zoom 
near the SAOs. The right part shows part of the trajectory in three-dimensional phase 
space and illustrates the early jumps of the stochastic sample path. 

Although a deterministic model is able to explain a some experiments (see, e.g., ex- 
perimental results by Hudson et al. |61]) it fails to accurately model realistic MMOs due 
to the presence of noise (see, e.g., the experiments by Dickson et al. [37]). Figure [l] also 
shows a sample path which is a stochastic version of the deterministic orbit perturbed by 
Gaussian white noise. Two important observations that can be made about the stochastic 
MMO pattern are 

• Part of the SAOs become indistinguishable from the random fluctuations so that 
counting SAOs below a certain amplitude is impossible. 

• The stochastic sample path typically jumps before the deterministic solution makes 
an LAO. 

In this work, we are going to provide rigorous formulations and the proofs of both obser- 
vations. 
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Our motivation (Ml) is to contribute to a better understanding of bifurcations in ran- 
dom dynamical systems by relating deterministic and stochastic methods for fast-slow 
systems. Within the last two decades substantial progress on deterministic fast-slow 
systems has been made. The analysis of hyperbolic fast dynamics has been completed 
in a series of works by Fenichel (see also [HB]) at the end of the 1970s; the theory 
focuses upon perturbations of normally hyperbolic critical manifolds to nearby slow man- 
ifolds. Near bifurcation points of the fast dynamics major developments in the 1980s 
used nonstandard analysis j30l HI] and asymptotic methods [l2l [HI [9l [86]. In the 1990s, 
two main geometric methods were introduced into multiple time scale systems. Switch- 
ing between hyperbolic fast and slow dynamics was analyzed using the exchange lemma 
[5 ^ I103j . The geometry of non-hyperbolic or singular dynamics has been uncovered using 
the blow-up technique [HI [10] . Since then the blow-up method has been applied success- 
fully for many singular-perturbation problems [HI UM [13] . In particular, folded nodes 
jlOU I107[ [76| [70] are a very interesting case since they occur already in generic fast-slow 
systems with one fast and two slow variables and have a highly nontrivial unfolding. Near 
a folded node the canard phenomenon occurs, i.e., orbits stay near a repelling manifold 
for a long time. This effect can generate MMOs [291 132]: see also the discussion for (M2) 
above and Figure [T] 

Noise acting on a system with multiple time scales can induce new phenomena such 
as early transitions \\M [M [HSl IHl [HI] and stochastic resonance [ISl [SD [HSl IH iHl HH] ■ 
The mathematical theory of fast-slow stochastic differential equations has mainly been 
developed during the last decade and is still far from being complete. Classical work on 
random perturbations of dynamical systems [H] , which mainly focused on large-deviation 
aspects, can be applied to situations with a time-scale separation exponentially large 
in the noise intensity [IHl UHl IS21 [SZl [SH]- A different approach, based on a detailed 
description of sample-path properties, applies to situations with time-scale separation and 
noise intensity of comparable magnitude [2n l22l [T9l 118] . This method led to a general 
theory for the behaviour of sample paths near normally hyperbolic invariant manifolds [231 
[21]. Other approaches include |;67J, which is based on moment estimates, and [96], which 
adopts the viewpoint of random dynamical systems in the sense of [6]. The associated 
methods and results have important applications in climate dynamics [HI [SU 1102^ [20] , 
the theory of critical transitions [951 E] , classical and quantum atomic physics [31 HI [2] 
and neuroscience [1051 [831 EH IHS [25] . In particular, canards in the stochastic FitzHugh- 
Nagumo system describing the action potential of neurons have been considered from 
the points of view of large deviations [361 EH] , and of convergence of sample-paths [HHj . 
Stochastic MMOs have also been considered in certain planar systems [HHl [HS] and in 
coupled oscillators [109] . 

The theory of stochastic differential equations with higher-dimensional singularities 
and multiple slow variables is not yet as advanced. Here we make a first step towards bridg- 
ing this gap between the generic higher-dimensional deterministic theory and stochastic 
sample-path analysis. The non-technical statements of our two main results are: 

(Rl) Sample paths near a folded node stay inside a tubular neighbourhood of an attract- 
ing deterministic solution. The neighbourhood is explicitly defined by the covari- 
ance matrix of a linearized process. The relation between the noise level, the time 
scale separation and a system parameter determines when small oscillations near a 
folded node become indistinguishable from noisy fluctuations. This relation can be 
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calculated explicitly to lowest asymptotic order. 



(R2) Sample paths near a folded node typically escape from a repelling deterministic 
solution earlier than their deterministic counterparts. The typical escape time can 
be determined rather precisely and depends on the same parameters as the relation 
in (Rl). The probability of observing atypical escape times can be shown to be 
small. 

Both results have important implications from theoretical and applied perspectives. 
In particular, we show how to control stochastic sample paths near a multi-dimensional 
bifurcation point. Therefore it is expected that the methods we develop have a much 
wider applicability beyond folded nodes, e.g. to singular Hopf bifurcations |52] or other 
stochastic bifurcation problems [SHE]. The precise quantitative estimates on the relation 
between noise level and a parameter controlling the number of SAOs are immediately 
useful in applications. Furthermore, the effect of early jumps could potentially regularize 
the complicated flow maps near a folded node [SH ES] and simplify the local-global de- 
composition of return maps [79] . 

The paper is organized as follows. In Section [2] we review the necessary theory for 
deterministic fast-slow systems and fix the notation. In Section [3] we state the known 
results about folded nodes and explain why they produce small-amplitude oscillations. 
In Section |4] we consider a variational equation around a special canard solution, called 
the weak canard. The solution of the variational equation can be transformed into a 
"canonical form" which allows us to prove a result on the spacing of canard solutions up 
cross-sections near or including the folded node. The proof is postponed to Appendix \K\ 
In Section |5] we develop the main setup for stochastic fast-slow SDEs and recall a result 
on attracting normally hyperbolic slow manifolds away from bifurcation points. Several 
notations that we use throughout are introduced as well. In Section [6] we give the rigorous 
formulation of our main results (R1)-(R2) for stochastic folded nodes. The result (Rl) 
on covariance tubes for the linearized process is proved in Appendix [Bj The influence of 
nonlinear contributions is dealt with in Appendix [C} The result (R2) on early jumps is 
proven in Appendix D Section [7] develops numerical simulations to visualize the analytical 



results. We conclude by giving a summary of parameter regimes and discussing the 
influence of early jumps on the global return mechanism and LAOs in Section |8| 
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2 Fast— Slow Systems 

We are only going to give a brief introduction to multiple time scale dynamics. A detailed 
reference covering many more topics is currently being written [80]; other, excellent ref- 
erences are [HH ED] for asymptotic methods and [ISj [TJ |66] for geometric methods. Many 
important discoveries were first made using nonstandard analysis [131 [38] : in particular. 
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many results we review in Section |3] were discovered by Benoit fTEj [TTl |T2] . However, we 
are not going to use any nonstandard methods and focus on the geometric viewpoint. 

A fast-slow system of ordinary differential equations (ODEs) is given by 



dy ■ / \ {^-^ 



where {x, y) G M™ x are phase-space coordinates, ^ & W are parameters and < e ^ 1 
represents the ratio of time scales. We shall assume that /, g are sufficiently smooth. By a 
rescaling we can change from the slow time s to the fast time t = s/e] this transforms 



— - x' = f{x,y,fi,e), , . 

dy _ „,/ _ ^^i^ „, ,, ^^ K'^-'^) 



\2.l\ to 

At = y = ^9[x,y,f^,e) . 

Remark: The more common notation for the slow time would be r but we shall reserve r 
for stopping times of stochastic processes; see Section [5] 

The first step to analyze fast-slow systems is to consider the singular limit e — )■ 0. 

(2.3) 



From (2.2) we obtain 

x' = /(x,?/,/i,0) , 
y' = 0. 

which is an ODE for the fast variables x where the slow variables y act as parameters. 



We call (2.3) the fast subsystem or layer equations; the associated flow is called the 



fast flow. Considering e — in (2.1) we find a differential-algebraic equation for the 
slow y-variables 

= f{x,y,fi,0) , ^2 4) 

y = g{x,y,iJ^,o) , 



called the slow subsystem or reduced system; the flow induced by (2.4) is called the 
slow flow. The slow subsystem is defined on the critical manifold 

Co:={(x,y)GM™+":/(x,y,/i,0) = 0}. 

Observe that Cq can also be interpreted as a manifold of equilibria for the fast subsystem. 
Note also that Cq does not have to be a manifold |73] but we only consider the manifold 
case in this paper. If the Jacobian matrix (D^/)(p) has maximal rank at p G Cq then the 
implicit-function theorem describes Co locally as a graph 

/io:M"^M'", f{ho{y),y,fi,0)=0 

near p. This allows us to write the slow subsystem more concisely as 

y = 9{ho{y),y,f^,o) . (2.5) 

We can strengthen the assumption on [Dxf){p) and require that it is a hyperbolic matrix, 
i.e., {D.j.f){p) has no eigenvalues with zero real part. In this case we say that Cq is 
normally hyperbolic at p. If all the eigenvalues of {D^f){p) have negative (positive) 
real parts we say that Cq is attracting (repelling) with respect to the fast variables. 
The following theorem shows that normal hyperbolicity is the key regularity assumption 
for fast-slow systems. 
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Theorem 2.1 (Fenichel's Theorem [13]). Suppose Mq is a compact normally hyperbolic 
submanifold (possibly with boundary) of the critical manifold Cq and that f,g & , 
1 < r < oo. Then for e > sufficiently small the following holds: 

(Fl) There exists a locally invariant manifold diffeomorphic to Mq. Local invariance 
means that can have boundaries through which trajectories enter or leave. 

(F2) Me has a H aus dor ff distance of 0{t) from Mq. 

(F3) The flow on converges to the slow flow as e — 0. 

(FJi^) Me is -smooth and can locally be given as a graph : — )■ M™". 

(F5) Me is normally hyperbolic with the same stability properties with respect to the fast 
variables as Mq. 

(F6) Me is usually not unique. In regions that remain at a fixed distance from the bound- 
ary of M^, all manifolds satisfying (F1)-(F5) lie at a H aus dor ff distance 0{e~^^^) 
from each other for some K > with K = 0{1). 

We call a perturbed manifold Me a slow manifold. Sometimes we refer to "the slow 
manifold" despite the non-uniqueness (F6) as it will be often irrelevant which of the 
(9(e~^/'^)-close manifolds we pick. 

A simple example where normal hyperbolicity fails is given by 

ex = y-x\ .2.6) 

y = jjL — X . 

The critical manifold Cq = {{x,y) E : y = x^} splits into three parts Cq = Cq U 
{(0,0)}UC5 where 

= Co n {x > 0} , and C^ = Cq n {x < 0} . 

Cq is attracting and Cq is repelling. At (x, y) = (0, 0) the critical manifold is not normally 
hyperbolic and has a generic fold singularity [71]; observe that (x, y) = (0, 0) is a fold 
(or saddle-node) bifurcation of the fast subsystem. Figure ^ illustrates the dynamics near 



the fold point of (2.6) 



To calculate the slow subsystem on Cq we could consider the two graphs x = hQ^y) = 
±y^ as suggested by (2.5). For (2.6) it is more convenient to differentiate y = x"^ implicitly 



with respect to s. This gives 

/i — X 



y = 2xx =^ X 



2x 



which shows that the slow flow is undefined at (0,0) if ^ 0. Fenichel's Theorem pro- 
vides slow manifolds and CJ. A major step in the theory of fast-slow systems was to 



consider the dependence of the dynamics of (2.6) on the value of n [H] El E3]. Note that 



for /i = the slow flow is well-defined and there is a special trajectory that passes from 
Cq to Cq and that a singular Hopf bifurcation [75] [27] occurs for = and < e ^ 1. 
The slow manifolds can be extended under the flow into the fold point region. Comparing 
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Figure 2: Planar fold near a singular Hopf bifurcation of (2.6); for (a) and (c) we have 
fixed e = 0.05. (a) fi = 0.1: The equilibrium point {xo,yo) = (/x,/^^) is determined as 
the intersection of Cq (grey) and the nuUcline {x = /i} (dashed black). The equilibrium 
is stable and the slow manifolds (red) and (blue) do not interact, (b) fi = 0: 
Only the slow flow on Cq (single arrow) and the fast flow (double arrows) are shown. Cq 
coincides with a maximal singular canard, (c) fi = —0.1: After the Hopf bifurcation the 
slow manifolds have "exchanged sides" suggesting an intersection for some /z near 0. 



Figure to Figure ^c) we expect that there is a parameter value /i for which the 
slow manifolds intersect/coincide. This intersection marks what has become known as a 
canard explosion [751 EH] . 



More generally, suppose that 7^ is a trajectory of a fast-slow system (2.1). Then we 
call 7e a maximal canard if it lies in the intersection of an attracting and a repelling slow 
manifold; for e = we also refer to 70 as a maximal singular canard. Canards in planar 
fast-slow systems are of codimension one whereas for higher-dimensional systems we do 
not need an additional parameter. In the next section we are going to focus on canards 
in three dimensions. 



3 Folded Nodes 

A general three-dimensional fast-slow system with one fast variable and two slow variables 
can be written as 

ex = f{x,y,z,i^,e), 

y = 9iix,y,z,iJ,,e) , (3.1) 
z = g2ix,y,z,ij,,e) . 



We assume that the critical manifold Cq = {{x,y,z) G M'^ : f{x,y,z,fj,,0) = 0} of (3.1) is 
a folded surface near the origin; suitable non-degeneracy conditions [321 llOlj are 

/(O, 0, 0, /i, 0) = , /,(0, 0, 0, /i, 0) = , 

/,(o,o,o,/i, 0)7^0, /,,(o,o,o,/i, 0)7^0, ^ ■ ' 

where subscripts denote partial derivatives. 
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Figure 3: 



Singular limit e = for the normal form (3.8) with /i 
red), the fold line L 



0.15. The attracting 

manifold Cq (red), the fold line L (green) and the repelling manifold Cq (blue) partition 
the critical manifold. Trajectories of the slow (fast) subsystem are indicated by single 
(double) arrows. 



The critical manifold again decomposes into three parts 

Co = U L U Co" , 

where Cq = CQr\{fx > 0} is repelling, Cq = Con{fx < 0} is attracting and L = Con{fx = 
0} is the curve of fold points; see Figure |3j Note that the assumption fy{Q,0,0, fi,0) ^ 



in (3.2) implies that the fold curve L can be locally parametrized by z. To obtain the 
slow subsystem we again differentiate /(x, y, z, /x, 0) = implicitly with respect to s 

xfx + yfy + zf^ = . 
This implies that the slow subsystem is 

-fy9i - fzg2 , 



z 



92 



(3.3) 



where all functions are evaluated for p = {x,y,z) G Cq and e = 0. On L the ODE 
(3.3) is singular but we can rescale time s t— )■ —s/f^ to obtain the desingularized slow 

(3.4) 



subsystem 



fy9l + fzg2 
-1x92 



Note that the time rescaling has reversed the orientation of trajectories of (3.3) on Cq 
and that (3.4) is a well-defined planar ODE. We define 

Kz)-= Uy9i + fz92)\p&L 
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and make the assumptions that 



m 

l{z) 



= 

7^ 0, 



Observe that (3.5) and /^^(0, 0, 0, /i, 0) 



for 2; 7^ 
imply that 



(3.5) 
(3.6) 



(x, z) = (0, 0) is an equihbrium 
point for (3.4) that hes on the fold curve. We say that {x,y,z) = (0,0,0) is a folded 
singularity. Points {x, y, z 0) ^ L are called jump points as trajectories have to make a 
transition from the slow to the fast flow at these points; the condition (3.6 ) is also called the 
normal switching condition. Generic folded singularities can be classified according to 
their equilibrium type into folded saddles, folded foci and folded nodes |10H [TT]. 
Folded nodes are the most interesting folded singularities. Without loss of generality we 
may assume that the folded node is stable for (3.4) with associated eigenvalues Ai, A2 for 
the linearization of (3.4). 



Proposition 3.1 ([29 | I107j ). Suppose that (3.1) satisfies (3.2), (3.5), (3.6) and that we 
are in the folded-node scenario. Then there exist a smooth coordinate change and a smooth 
change of time which bring (3.1) near {x,y,z) = (0,0,0) into the form 



ex 

y 

z 



y — x"^ + 0{yx^, x^, xyz) + eO{x, y, z, e) 
-{fx + l)x-z + 0{y,e,{x + zy) , 

2 ' 



(3.7) 



with Ai = —fi and A2 = — 1 being the eigenvalues for the linearization of (3.4). 



We shall show in Section |4| that the terms in (3.7) denoted by 0{-) are indeed higher- 
order for the analysis near the folded node. Hence we can work with the normal form 



ex 

y 

z 



y~x\ 
-(/i + l)x 



(3.8) 



The critical manifold of (3.8) is 

Co 

It splits into three components 



{{x,y,z) eR^ ■.y = x'} . 



Co 



where Cq = Co H {x > 0} is attracting, Cq = Cq fl {x < 0} is repelling, and L is now a line 
of fold points. Parametrizing over the slow variables, we can also write Cq'^ = {x = ±^/y}, 
i.e., h^'^^iy^z) = 

Again differentiating y = x^ implicitely, we see that the slow subsystem is given by 



2xx 

z 



-(/i + l)x — z 



2 ' 



(3.9) 



and the desingularized slow subsystem (see also Figure |4]^a)) is of the form 



X = — (/i + l)a; — z 
z = fix . 



(3.10) 
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Figure 4: (a) The desingularized slow flow (3.10) is sketched for some /i > with a stable 
node at the origin, (b) The slow flow (3.9) is illustrated. The strong eigendirection (solid 
red) deflnes the strong singular canard 7q; the weak eigendirection (dashed red) deflnes 
the weak singular canard 7^. 



The system (3.10) is linear with an equilibrium point at (x, z) = (0, 0), and the eigenvalues 
are (Ai, A2) := (—1, — /u). We assume from now on that 

/iG (0,1) 



so that (0, 0) is a stable node for the desingularized slow subsystem (3.10). Hence we also 
denote the eigenvalues as 



Ai = -1 



A. 



and 



A2 = -/i =: A^ 



to emphasize the strong and weak eigendirections. Note that fi = A^/Ag precisely rep- 
resents the ratio of eigenvalues and attains all resonances fi~^ G N for G (0,1). The 
associated (unnormalized) eigenvectors are 



and 



(3.11) 



which also represent directions for two maximal singular canards; see Figure |4} Observe 
that the singular strong canard 7q and L bound a funnel region on Cq of trajectories 
that all flow into the folded node; see Figure |4} The funnel region has an opening angle 
cos^^{jj,/ a/1 + /i^) which converges to 7r/2 as — ?■ 0. The funnel on Cq is located in the 
{x > 0, z < 0}-quadrant and the singular weak canard 7,^ is given by the anti-diagonal 
{z = —x}. 



Theorem 3.2 ( |101l fTTl [T5 | 1107] ). Suppose (3.1) has a generic folded node (i.e. Proposi 



tion 3.1 applies). Then for e > sufficiently small the following holds: 



(CI) The singular strong canard 7q always perturbs to a maximal canard 7^. // /i^^ N, 
then the singular weak canard 7^ also perturbs to a maximal canard 7^. We call 7^ 
and 7^ primary canards. 
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(C2) Suppose k > is an integer such that 



2k + l < fx'^ <2k + 3 and ^"^ ^ 2{k + 1) . 

Then, in addition to ■yp'^ , there are k other maximal canards, which we call secondary 
canards. 

(C3) The secondary canards converge to the strong primary canard as e — > 0. 

(C4) The primary weak canard of a folded node undergoes a transcritical bifurcation for 
odd G N and a pitchfork bifurcation for even fi~^ G N. 

We emphasize the results (C1)-(C3) which will be of major importance for our stochas- 
tic analysis; (C4) describes the behaviour near resonances and will not be considered here. 
The next theorem provides a geometric viewpoint for the generation of maximal canards 
near a folded node. We say that a twist corresponds to a half rotation (i.e. a rotation 
by an angle of vr). 

Theorem 3.3 f pUTl [TUT] ). Assume 2k + 1 < jj,'^ < 2k + 3, for some k e N, and 
H'^ 7^ 2{k + 1). Then the following holds: 

(C5) The primary strong canard 7o twists once around the primary weak canard 7^. 

(C6) The j-th secondary canard 7^, 1 ^ j ^ k, twists 2j + 1 times around the primary 
weak canard 7^ . 



(C7) The twisting/rotation occurs in an 0{^/e) neighborhood of the folded node for (3.1). 



In particular, the slow manifolds and start to spiral near the folded node creat- 
ing transversal intersections away from resonances. For visualizations of these manifolds 
in several different contexts see |33l [32l 13] 



In Figure |5j we show the slow manifolds for (3.8) near the folded node on the cross- 
section {z = 0}. The manifolds have been computed by forward integration and using 
the symmetry 

(x, y, z, s) ^ (-X, y, -z, -s) . (3.12) 

The center of rotation is the weak canard 7,^. Since /i = 0.08 we know by (C2) that there 
are five secondary canards in Figure [5j Five intersections are indeed detected numerically 
but 7q'^ are very close to 7,^ on {z = 0}. All secondary canards 7^ approach 7^ when 
z < near the folded-node region on . The canards move away from each other for 
z > (see Figure [3]). The next theorem shows that the maximal canards organize the 



rotational properties of trajectories passing through a folded-node region. 
Theorem 3.4 ([29j). Fix two sections 

Si := {(x, y,z) eR^ -.y = Ki} for some < Ki = 0{1) 

:= {(x, y,z)eM? ■.y = K2t} for some < K2 = C(l) 



for the system (3.8). Consider the intersection points of maximal canards in S-'^flC". Let 
k be the number of secondary canards. Then, for e > sufficiently small, the following 
holds: 



11 




Figure 5: Canards and slow manifolds near a folded node for (3.8) with (/i, e) = 
(0.08,0.01) on the cross-section {z = 0}. The primary canard 7^ and the first two sec- 
ondary canards 7^'^ are labeled. We also show a zoom near the primary weak canard 7,^. 
All maximal canards are indicated by black dots. 



(C8) The secondary canards are 0{e^^ ^^^'^) close to the primary strong canard. 

(C9) There exist {k + 1) so-called sectors of rotation Ij, 1 ^ j ^ A; + 1, between the 
two primary canards labeled in increasing order starting from the strong primary 
canard. The size of the sectors Ij for 1 ^ j ^ k is 0{e^^~^^^'^) while the size of the 
sector Jfc+i is 0{1). 

(CIO) The Poincare map from T,^ to is a contraction with rate 0{e^^^'^^^^'^'^^) . 

(Cll) All maximal canards are separated by 0{^/e) in their z-coordinate on T? . 

Remark: We note that the results from Theorems |3.2 3.3 and 3^ also extend to higher- 
dimensional fast-slow systems with at least two slow variables and at least one fast variable 
|108[ |29] but that they do not provide a detailed analysis of canards beyond the section S^. 



Theorem 3.4 provides sectors of rotation that organize the twisting of trajectories near 
the folded node. Once we know which sector an orbit enters we can predict the number of 
oscillations. Note that the oscillations can be classified as "small oscillations" due to (C7). 
Global returns can induce so-called mixed-mode oscillations (MMOs) which are found 
in a wide variety of applications; see [32] for a review of MMO mechanisms in multiple 
time scale systems. 
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4 Canard Spacing 



Theorem 3^ describes the spacing of maximal canards away from the folded-node region. 
Since we are also interested in their spacing on the cross-section {z = 0} depending on 
H we need a refined analysis near the folded node. The key component in the proofs of 
Theorems 3.2, 3.3 and 3.4 is a rescaling of (3.8) near the folded node 

(x, y, z, s) = {y/ex, ey, y/ez, y/es) (4.1) 

which can also be interpreted as a blow-up transformation |10H 11071 129] . We shall not 
introduce the blow-up method here but restrict ourselves to the analysis of the rescaled 
system 

X = y-x^ + O(v^) , 

V = -{fi + l)x-z + 0{^e) , (4.2) 

I = 



Therefore the (9(-)-terms in (3.7) are indeed of higher order for the analysis near the 
folded node. Neglecting the small e-dependent terms and dropping the overbars in (4.2) 
for notational convenience yields 



X 

y 



y — X 

2 ■ 



l)x 



(4.3) 



The ODEs (4.3) are our main focus of study in this section. Note that (Cll) in Theo- 
rem 3.4 implies that the maximal canards are all (9(l)-finitely separated for (4.3) when 



they are an 0(l)-distance away from {x,y,z) = (0,0,0). Observe that we can always 
solve the last equation 

= ^{s - So) + Zo 

where zo = z{so) denotes the initial 2;-coordinate of the trajectory at the initial time 
s = So- Hence we can view 2; as a time variable and re-write (4.3) if necessary as a planar 
non-autonomous ODE 

2y - 2x^ , 
-2{fi+l)x - 2z . 



(4.4) 



Our first goal is to quantify the intersections of canard solutions with the section {z = 0}. 
Therefore we are going to focus on the analysis of orbits arising as perturbations of slow 
subsystem trajectories inside the funnel on Cq and assume 

20 ^ ^ ^ . 



Let {x* (z) , y* (z)) be any solution of (4.4) and set 



u 



to derive the variational equation 



[X 



X ,y-y 



-4mix* + 2u2 - 2ul 
-2(/x + l)ui . 



(4.5) 



A key observation by Benoit [13] was that there are some special solutions to (4.3). 
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Lemma 4.1 (^13j). The ODE (4.3) admits two polynomial solutions 

{x{s),y{s),z{s)) = Q^'^^^ + ^'f* 
with X G {— /U, —1} corresponding to the two primary singular canards. 



(4.6) 



Lemma 4.1| can be checked by direct differentiation of (4.6). We know that the weak 
canard is the center of rotation, and from (4.6) with A = — /i we find that the variational 
equation around the weak canard is given by 

du 



dz 



Az 2 
-2(^ + 1) 



Ui 
U2 



+ 



-2ul 




(4.7) 



We are interested in the detailed interaction of other maximal canards witli the weak 
canard. 

Proposition 4.2. If fi > is sufficiently small and z < is bounded away from then 
solutions u = u{z) of (4.7) are attracted exponentially fast to {ui = = U2}- 



Proof. Augmenting the variational equation (4.7) by i = gives an autonomous fast-slow 
system with two fast variables u = {ui,U2) and one slow variable z (since fi is sufficiently 
small). The critical manifold is 

{{Ui, U2, Z) ■.Ui= 0, U2 = 0} . 

Computing the linearization with respect to the fast variables for yU = gives a matrix 

' iz 2 
-2 

with eigenvalues 2{z± ^ z'^ — 1). Hence if 2; < both eigenvalues have negative real parts 
and the result follows from Fenichel's Theorem. □ 



For < /i ^ 1, Proposition 4.2 allows us to reduce the study of the nonlinear varia- 
tional equation (4.7) to a linear one by dropping the higher-order term —2u\. This yields 



the linear non-autonomous ODE 

du 



[I 



dz 



4z 2 
-2(/i + l) 



u = A{z)u . 



(4i 



--■.A{z) 



In particular, we must show what happens to solutions u = m(ui, M2) near the weak canard 
and when z is not bounded away from 0. 



Remark: The variational equation (4.8) has been analyzed llOlj by re-writing it as a 
second-order equation 

d^ui dui 



ds2 



— jJS- 



ds 



+ ui=0 



Benoit [13] observed that using the time rescaling s = s/^/JI of (4.9) one gets 





d^Mi _dMi 1 
dp ds 



(4.9) 



(4.10) 
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which is referred to as Weber equation or Ricatti-Hermite equation. The ODE (4.10) 



has exphcit solutions in terms of Hermite polynomials (see [Ij, p. 781). Then the asymp- 
totic properties of Hermite polynomials can be used to draw conclusions about the exis- 
tence of maximal canards. 



We develop two alternative ways to describe the variational equation (4.8) directly. 



This provides new quantitative information about canard solutions and also gives esti- 
mates for the spacing of canards on cross-sections near the folded node. Our first approach 
in Section [ITT] uses averaging and provides a formal result. This result is then stated in 



Section 4^ and proven in Appendix |X] using coordinate changes that are motivated by 



the formal calculation. 



4.1 Averaging 



As shown in (4.9) we can consider the variational equation (4.8) as the second-order 
equation 

■ill + Ui = fi'Siii . (4-11) 



The form (4.11) suggests to view the problem as a damped oscillator; we assume that 

2 



fx 



< So ^ s < . 



(4.12) 



Note that /is = —2 is precisely critical damping and for fj,s e (—2, 0) the oscillator 
is underdamped. Therefore we expect that (4.11) describes harmonic oscillations with 
damping/contraction that is non-uniform in time. We make a change to polar coordinates 

(mi, ill) = (r(s) cos(s + ipi^s)), — r(s) sin(s + ip^s))) 

which yields the ODEs 



r = fisr sin^(s + ip) , 



To simplify the analysis, consider the time rescaling s 



(4.13) 



s. This converts (4.13) to 



dr 

ds 



-f sin2(-v^ + ^) 



cos(— y— s + ip) sin(— V— s + ip) . 



(4.14) 



We consider (4.14) on each time subinterval 

/, := [-((j + 2)7r)2, -{j7r)% for j G {0, 2, 4, 6, . . .} 



by viewing (4.14) as a vector field on IR+ x (Ij/ ~) where the equivalence relation ~ 
identifies the endpoints of Ij. Then the vector field is in the form for averaging |106j . The 
averaged equations are 







ds 








ds 


2\I, 



sin^(— -|- ijjj)ds 



(2j + 2)7r + sin(2?/'j) 



2(4j + 4)7r 



cos 



-s + ipj) sin(— -\/^ + il>j)ds 



fITf cos(2^/'j) 
"2(4j +4)7r2 



(4.15) 
(4.16) 
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In particular, we find that if we take a formal limit j — oo the equation for the radius is 



ds 



^0 

4 



(4.17) 



Remark: Observe that one could also view (4.15)-(4.16) as an autonomous vector field 
and formally average over the angle ipj to get 



U + 1) 
ds (4j + 4) 

and then take the limit j oo. 



ds 







The solution of (4.17) is given by 



roo(5) = r„o(s~o)e-i'^(^'-^''') =r„o(so)ei'^(^'-^o) =r„o((2//x)^o)e(^'-^o)/^ (4.18) 
which shows that the leading-order behaviour of the solutions to the variational equation 



(4.8) for 2; < consists of a contraction towards the weak canard, given by (4.18), com- 



bined with a rotation. If we assume that ipo{0) = 0, i.e., the rotation ends at angle on 



section {z = 0}, then solving (4.16) yields 



tpo{s) = —tan ^ 



tanh 



;i+j> 



In principle we can now calculate ^/'o(— (2vr)^), use this result as an initial condition for ipi 
and then repeat the process to get a very detailed description of the rotational properties 
of trajectories near a folded node. 

4.2 Diagonalization 

To give rigorous arguments instead of the above formal calculation, we start by considering 



the variational equation in first-order form (4.8). The matrix A{z) has eigenvalues 

2z ± 2iu{z) , where ijj{z) = ^/l — z"^ + fi . 
We assume that |z| < 1 so that uj{z) is real and bounded away from 0. Furthermore A{z) 



has trace Az. We expect that the solution u{z) for (4.8) consists of a contraction and a 
rotation for 

1< zo^ z <0 . (4.19) 



Observe that (4.19) corresponds to the condition (4.12) 



Theorem 4.3 (Canonical form). There exists a matrix 

1 



Siz) 



—z + uj{z) —z — u{z) 



(4.20) 



such that the coordinate change u{z) = S{z)u{z) transforms the variational equation (4.8) 
into canonical form 

du ( a(z) w{z, , 2ij 



^ dz \-w(z) a{z) ^ 
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where 



a{z) = 2z + 0{n^) , 
w{z) = 2uj{z) + 0{fi) . 



As a consequence, we can write the solution of (4.8) in the form 

u{z) = e""^'''"^/" S{z)U{z, zo)S{zo)'^u{zo) 
for \z\ < 1 and \zo\ < 1, where 



(4.22) 
(4.23) 

(4.24) 



a{z, zq) = / a{s) ds = z'^ — Zq + 0{iJ?) 

J Zo 

and U{z, Zq) is the orthogonal matrix 

TKt 7\- ( cos{^{z,s)/fi) sm{ip{z,s)/^I)\ , . _ C , . A.^niu\ 



(4.25) 



The proof is given in Appendix |Aj We remark that Theorem 4^ easily extends to the 
hnearization around solutions other than the weak canard. 

4.3 The Distance Estimate 



From the solution of the variational equation in Theorem 4^ we can give asymptotic 
estimates on the distance of the secondary canards to the weak canard on the section 

SO := {{x,y,z) G : ^ = 0} . 

Theorem 4.4 (Canard spacing). The distance of the k-th secondary canard 'y^ to the 
weak canard 7"" on S° is given by 0{e~^^'^''~^^'^ ^) as fi where < cq = 0{1) is a 
constant. 



Proof. By Theorem 3.3, part (C6), the k-th. secondary canard makes {2k + l)/2 twists 



around the weak canard for z going from zq to 0. Using this fact together with the 



rotational part (4.25) of the solution (4.23) in Theorem 4.3 we find the condition 

ip{0,zo) 



/i 







w(s) ds 



2k + 1 



TT- 



(4.26) 



to leading order. Using concavity of the integral z H- u{s) ds for z G (—1, 0), we find 



TT 



1 + /i ^ / w(s) ds ^ — Zqa/I + 

J Zo 



(4.27) 



which enables us to approximate the integral in (4.26). Hence as /i — )■ we get that 

TT 



Zq G 



i2k + l)fi,--{2k + l)fi\ =: [4'\4'^] 



(4.28) 



where zj^'^ is the starting-point estimate for maximal canards that start to spiral around 
7^" near 2; ^ — 1 and z^^ the estimate for maximal canards that start to spiral around 7"' 
near z '^0. In particular, we find that zq = 0{{2k + l)/i). Note that the contraction term 



towards 7^" in (4.23) for —1 < ^ -2 < has order 



The result follows upon evaluating (4.29) on S and substituting (4.28) 



(4.29) 

□ 
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5 Stochastic Fast-Slow Systems 



As for deterministic fast-slow systems we only give a brief introduction to stochastic 
fast-slow systems. For a detailed introduction to stochastic multiple time-scale dynamics 
consider EO]; we are also going to assume standard results on stochastic differential 
equations (SDEs) [92l EH] • AH SDEs and stochastic integrals in this paper are considered 
in the Ito interpretation. 



We associate to the deterministic fast-slow system (2.1) a stochastic fast-slow 
system given by 

dxs = ^f{x,y,fi,e)ds + -^^F{xs,ys)dWs , 
dys = g{x,y,fi,e)ds + a'G{xs,ys)dW, , 

where F : M'^+" G : {Wsjs^o is a fc-dimensional standard Brow- 

nian motion on some probability space (fi, J-", P) and a,a' > are parameters controlling 
the noise level 

We also define p := a' /a and assume that p is bounded above and below by positive 
constants. The initial conditions (2:0,2/0) := {xsQ,yso) are chosen to be square-integrable 
and independent of {Wsj^^sg. Furthermore we assume that f,g,F,G are sufficiently 
smooth and satisfy the standard assumptions that guarantee the existence and pathwise 



uniqueness of strong solutions {xs,ys) to the SDK (5.1). Note that these conditions also 
imply the existence of a continuous version of {xs,ys). The law of the process {xs,ys), 
starting at time sq in (xo,?/o)) is denoted by P^o,(a;o,?/o) a^^j ^^j^g corresponding expectation 



by E,^o,i^o,yo) _ q^j- approach to understand the dynamics of (5.1) is to analyze the time 
a sample path spends in a given Borel-measurable set A C M'"+". Suppose (xo,?/o) £ A 
and define the first-exit time from A as 

r4 := inf{s G [so,oo) : {xs,ys) ^ A} . 

In this paper we are always going to choose sets A so that r4 is a stopping time with 
respect to the filtration generated by {{xs,ys)}s^so- 

Our first goal is to state an analog of Fenichel's Theorem. This theorem is going 
to describe the typical spreading of sample paths near an attracting critical manifold 



bounded away from folded singularities. Suppose the deterministic version of (5.1) with 
a = = a' has a compact attracting normally hyperbolic critical manifold 

Go = {ix,y) G M"+" : x = hoiy),y E Vo} . 



Let Ce be the slow manifold obtained from Fenichel's Theorem 2A_ Our strategy is to 
construct a neighborhood B{r) for G^ that contains the sample paths with high probability 
[2SII21]. Define 

^s = Xs- K{ys) ■ (5.2) 
Observe that measures the deviation of the fast components from G^. Applying Ito's 
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formula to (5.2) gives: 



de. = dXs-Dyh,(y,)dys + 0{{a'Y)ds (5.3) 
= - [fiheiVs) + 6, Vs, fJ-, e) - ^Dyhe{ys)9{K{ys) + 6, Vs, fJ', e) + C(e(o-')^)] ds 



ye 

From now on, we suppress the arguments fi and e for brevity. Consider the hnear approx- 
imation of (5.2) in ^s, neglect the Ito term 0{e{a')^) and replace yg by its deterministic 



version yf'^* to obtain 



(5.4) 



where the two matrices and F° are defined as 



Aiy) = DJ{K{y),y)-tDyK{y) D,g{K{y),y) 

F^iy) = F{K{y),y)-pV~^DyK{y)G{K{y),y) ■ 



Observe that Ao{y) = Dxf{ho{y),y) and F^^y) = F{ho{y),y). To solve (5.4), we pick an 
initial condition on the slow manifold (^0 5?/o°*) — (OjZ/o'^*)- Then the solution of (5.4) is 
the Ito integral 



.0 
Ss 



cr 



f/(.,r)Ff(yf*)diy, 



where U{s, r) denotes the principal solution of the homogeneous linear system ez> = 
A^{yf'^)v. If we fix a time s then is a Gaussian random variable of mean zero and 
covariance matrix 

Cov(C°) = - fu{s,r)F',{yf')Fe{yTfU{s,TfdT . 
£ Jo 

Note carefully that Xg := o"^^ Cov(^^) does satisfy a fast-slow ODE given by 



eX = A,{y)X + XA,{yY + F^{y)F^{yY , 

y = 9{he{y),y) ■ 

The system (5.5) has a critical manifold 5*0 given by the equation 

Ao{y)X + XAoiyf + F^{y)Fl^{yf = . 
By the remarks above we see that this is equivalent to solving 

{DJ)ihoiy), y)X + X[{DJ){h,{y),y)]^ + F{ho{y).y)F{ho{y), yf 
Again this manifold can be locally described as a graph 

Sl = {{X,y)eW-'--:X = H,{y)} . 
The next lemma states that Sq is normally hyperbolic and attracting. 



(5.5) 



. 



(5.6) 
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Lemma 5.1 ([ID])- Let Mi and M2 be square matrices of dimension m with eigenvalues 

. . . Ai^m and X2,i, ■ ■ ■ ^2,m, respectively. Then the linear map L : 
defined by 

L{X) = MiX + XM2 
has m? eigenvalues given by {Xi^i + \2,j} for i,j G {1,2,..., m}. 



Therefore Fenichel's Theorem |2 . 1 1 provides us with a slow manifold 

Si = {iX,y) eve : X = H,{y) = H,{y) + 0{e)} , 

which we will now use to describe the typical spreading of sample paths. 

Theorem 5.2 ([SS','^^). Suppose the norms \\H^{y)\\ and \\H'[^{y)\\ are uniformly bounded. 
Define the neighborhood B{r) around the deterministic slow manifold as 

B{r) := {{x,y) eV:{[x- K{y)], H:\y)[x - K{y)]) < r^} . 

Then, for e > and a > sufficiently small, sample paths starting on remain in B{r) 
with high probability; in particular, we have 

p5o,(xo,yo)|^^^^^ < s A ri,J < e) q-^^"^'/^"' . (5.7) 

with K12 > 0. 

Detailed discussions of the factors Ki{s, e) and K2 can be found in [231 [2l]. partic- 
ular, Ki{s, e) grows at most like s^, while K2 does not depend on time and can be taken 
close to 1. The probability in (5.7) thus remains small on long time spans as soon as we 
choose r ^ cr. This implies that sample paths stay for exponentially long times near an 
attracting slow manifold before they jump away unless they come close to the boundary 
of Co before, i.e., the ^/-coordinates leave the set Vq. 



6 Stochastic Folded Nodes 



6.1 B low-Up 



Proposition 3.1 shows that Equation (3.8) is a normal form for deterministic fast-slow 



systems with a folded node. We study the associated SDE 



dxs 



(1) 



[-(/i + l)xs - Zs] ds + a' dW, 



(2) 



(6.1) 



fd. 



where wj^\ Ws^' are independent standard Brownian motions; to simplify the notation 
we also define Ws := (W^.W^'^^f. Since z plays the role of a time variable we do not 
add noise to the z-component. We will always assume that the noise terms are of equal 
order, i.e., that p = a' jo is bounded above and below by positive constants. Note that 



K2) 



(6.1) fits into the framework of a general fast-slow SDE (5.1) with 



F{x,y) = {\ 0) and G{x,y) = {^ l) 
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We apply the rescaling/blowup given by (4.1) to (6.1) to get 



/2j , 



e [-(/i + l)a;, - z^] ds + a' dWlt/2, , 



K2) 



We use the scahng law of Brownian motion, divide the first equation by y/e and the second 
one by e and drop the overbars for notational convenience, obtaining 



dxs 

dys 

dz. 



(y,-a;2)ds + ^diyP, 

[-{^i + l)x,-z,] ds+^.dW^''^ 

fd.. 



(6.2) 



Therefore rescaling the noise intensities as 



cr, O" 



removes e from the equations and yields 



dxs 

dys 
dz^ 



(1) 



(6.3) 



{ys-xl)ds + (jdW, 
[-(/i + l)xs - z,] ds + a' dW, 



(2) 



(6.4) 



fd.. 



where the overbars from (6.3) have again been dropped. To study (6.4) we consider the 



variational equation around a deterministic solution (which may be the weak primary 
canard, a secondary canard, or even any other solution with initial condition (a;o, yo? -^o) 
sufficiently close to Cq ). This approach is a generalization of Section where we considered 



the variation around the slow manifold. Viewing (6.4) as a planar non-autonomous system 
with time z = (/i/2)s and setting 



{^zi Hz 



we get the SDE 



diz 
dr]z 



H'n^-ez-'^xf'Qdz + a 



-l{^^ + l)i.dz + a'^^ 



^dW, 



(1) 



dW, 



(2) 



(6.5) 



6.2 Covariance Tubes 



We linearize (6.5) and denote by = (^°,?7°)^ the solution which satisfies the SDE 



dC 



1 / -4xf* 2Vo , , ^ /^v^ , ,,,, 



=:A(a:dct) 



,F0 



Then (6.6) is solved by the following Gaussian process 



C° = Uiz, zo)C + ^ / Uiz, r)F' dWr 



(6.6) 
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where U{z, r) is the principal solution to the deterministic homogeneous non-autonomous 
hnear system fj,( = A{x'^'^^)(. The two-by-two covariance matrix Cov(C°) =: Cov{z) is 
given by 

2 pz 

Cov{z) = — U{z,r){F^){FYU{z,rfdr. (6.7) 

Jzo 

Differentiating V{z) := cr^^ Cov(z) we find that it satisfies the ODE 

dV 



/i 



A{xf')V + VA{xfY + (F°)(F°)^ 



dz 



{61 



with initial condition V{zo) = 0. The following result describes the behaviour of the 



solutions of (6.8), in particular as z approaches 0. 



Theorem 6.1 (Behaviour of the covariance matrix). Fix an initial time Zq < 0. There 



exist constants c+ > c_ > such that, for a// < ^ 1, the solution of (6.8) with initial 
condition V{zq) = satisfies 



and 



^ ^ V^ll(^),V^22(^) ^ ^ 
^^V,,iz),V22iz)^^ 



|V^12(^)| = \V2l{z)\^C+ 
\V22{z)-Vu{z)\^C+ 



for zo + 0{fi\\ogfi\) ^z^ , (6.9) 

for-^^z^^, (6.10) 



for Zq ^ Z ^ ^/JI , 



for Zq ^ Z ^ 



//i 



(6.11) 



Furthermore, let V{z) be any solution of ( |6.8 ) with positive definite initial condition 
V{zo). More precisely, we require both V{zo) and V{zo)~^ to have all elements uniformly 



bounded in /z. Then the matrix elements ofV{z) satisfy (6.10) and (6.11), and (6.9) holds 
for all z G [zq, —y/JI]. 



The proof is given in Appendix [B| where we also give some additional information on 
how the covariance can be approximated by asymptotic expansions. 



Equations (6.9) and (6.10) show that the variances of and rj^ grow like o''^/\z\ up to 



time — i/yU, and then stay of order cr^/ y/Jl up to time ^/JI. Thus we expect the fluctuations 
of stochastic sample paths around the deterministic solution to increase when the fold at 
z = is approached. The restriction z ^ zq + 0{^\\ogfi\) is due to the fact that the 
variances are initially equal to zero, and need some time to build up. 



Equations (6.11) show that the covariance of and r]^ remains bounded, of order o"^, 
up to time 



//i, and the same holds true for the difference between the variances. This 
implies that fluctuations become more isotropic as z approaches 0. 



We now turn to the analysis of the full nonlinear SDE (6.5) satisfied by the difference 
(z = {^z,Vz) between stochastic sample paths and deterministic solutions. This SDE can 
be written in vectorial form as 



dc. = - [A{xi^')C + KCz 



a 



dz + —F^ dW, 



(6.12) 
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where A and F have been defined in (6.6), and &(C) = (— , 0) denotes the nonhnear 



term. We expect the covariance matrix of Q to be close to Cov(^°) = a'^V{z), where V{z) 



is the solution of (6.8) with initial condition V{zo) = 0. Thus sample paths should be 
concentrated in a tube surrounding the deterministic solution, with elliptical cross-section 
determined by V{z). Note that the elliptical cross-section becomes close to circular as z 
approaches 0, since the variances Vu, V22 are then of larger order than the covariance V12 
and the difference V22 — Vu. 

The fact that V{zq) is not invertible causes some technical complications. Therefore, 



in the following we let V{z) be the solution of (6.8) with an initial condition V{zo) 
which is positive definite. Observe that the difference between V{z) and V{z) decreases 
exponentially fast. We fix a 2:0 < and define the covariance-tube as 

B{r) = {{x,y,z): z^ ^ z ^ {[{x,y) - {xi^\ yt^%V (z)-' [{x , y) - {xf\yT')]) < r'} . 

(6.13) 

The cross-section of B{r) at any plane {z = const} is an ellipsoid whose axes are deter- 
mined by V{z), while the scaling parameter r controls the size of the tube. 

Theorem 6.2 (Concentration of sample paths in the covariance tube). There exist con- 
stants Ao, ro, /io > such that for all < A < Aq, all a < r < r^ji^/^ and all /i < /iq, 

P{t-bm < ^} ^ C^{z, zo) e-^'^'/'-' (6.14) 

holds for all z ^ ^/JI, where the exponent kq satisfies 

Ko = l- C(A) - C(r/i"3/^) - 0{ayr^) , (6.15) 

and the prefactor is given by 



const / r ^ ^ 
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C+{z,Zo) = l + =\^-] I xf'ds. (6.16) 



The proof is given in Appendix [C} This result shows that the probability that sam- 
ple paths leave the covariance-tube before time z is small, provided we choose r ^ 
a{\og{C+{z, zo))Y^'^. Indeed, for these values of r, the exponential term dominates the 
polynomial prefactor. We thus say that sample paths are concentrated in the covariance 
tube B{r) for r slightly larger than a, or that the typical spreading of sample paths is 
given by B{a). 

The condition a < r < tq^^/'^ implies that the theorem only applies to noise intensities 
smaller than (9(/i^/'^). What happens for a ^ /i^/^ is that fiuctuations become large 
already some time before the fold line is reached, which completely smears out the small 
oscillations present in the deterministic case. In fact, it is possible to show that if a ^ Z^^''^, 
the bound (6.14) still holds true for r ^ ro|-2p''^, and thus sample paths are localised up 



to times z <^ — cr^/^ 



6.3 Small- Amplitude Oscillations and Noise 

Our knowledge of the size of the covariance tubes now allows to determine when they start 
to overlap; from this we can deduce consequences for the existence of small-amplitude 
oscillations (SAOs) in the presence of noise. Recall from Theorem 3.2 that there are two 
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primary canards 7^'"^ and K — 1 secondary canards 7^ for 2K — 1 < fi^^ < 2K + 1. 
We want to consider the rotations around 7^ and denote the strong canard by 7°, i.e., 
k = 0. By Theorem 4.4 the distance from the weak canard is given by (9(exp(— co(2A; + 
for some positive constant Cq G [7r/4, 1]. Theorem 6.1 imphes that on z = the 
width of the covariance tubes is given by 0{jj^^^^a). Therefore the covariance tubes 
for sample paths starting on the A;-th canard start to overlap with the weak canard for 
a ~ /i^/'' exp(— co(2A; + We define the functions 



^l/4g-co(2fc+l)V_ 




Figure 6: (yU, cr)-parameter plane, (a) The curves crfc(/i) are shown for cq = 1. Above the 
k-th curve we cannot distinguish the oscillation induced by the fc-th canard from noisy 
fluctuations. The regions of different numbers of SAOs are delimited by the curves ak{fi) 
and the existence requirements for the k-th canard. The number of visible canards is 
indicated in rectangular boxes, (b) A zoom of (a) is shown that illustrates the structure 
of thin regions as /i — )■ and a 0. 

The previous considerations yield the following result. 

Corollary 6.3 (Noisy SAOs). On the section {z = 0} the covariance tubes of the k-th 
canard overlap if 

(y > cxM . (6.17) 

Therefore, depending on the noise level a, the canard number and the parameter n, deter- 
ministic SAOs with 2fc+l twists become indistinguishable from noisy fluctuations whenever 



(6.17) holds. 



In Figure |6] we show the curves crfc(/i) and indicate in which regions of the (/i, cr)- 
parameter plane one can distinguish which number of canards. The curves (Tfe(/i) and the 
existence conditions for canards enclose bounded regions where precisely + 1 canards 
can be distinguished which yields {2k + l)/2 twists up to the section {z = 0}; see also 



Theorem 3.3 We can also study Figure p\ for fixed a. In this case decreasing fi first 



increases the number of visible SAOs and then decreases it again. 
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6.4 Early Jumps 



We now turn to the behaviour for times z > ^/JI. For definiteness, we let (x^'^*, y^^*) = 
{—z, z^ — /i/2) be the weak canard solution, and define the set 

V[r^) = {{x,y,z):z^ ^, {x - xf'f + {y - yf'f < v'z] . (6.18) 

'D{ri) is a tube centred in the weak canard, whose width grows like ^/z. The following 
result, which is proved in Appendix |Dj shows that sample paths are unlikely to stay for 
very long in T>{'r]). 

Theorem 6.4 (Escape of sample paths from the primary canard). There exist constants 
K = K{ri) > 0, Co > and 71,72 > such that, whenever a\\oga['^ ^ A^^''^, 

^{rv{^) > z] ^ Co|log(T|^2e-''(^'-'^)/(^|i°s-l) . (g.19) 
The probability that a sample paths stays in 'D{ri) thus becomes small as soon as 

z > V/i|log(T|/K . (6.20) 

Unless the noise intensity a is exponentially small in /i, the typical time at which sample 
paths jump away from the canard is slightly (that is, logarithmically) larger than ^JJi. 



7 Numerics and Visualization 

In this section we briefiy discuss how to compute canard solutions and their associated 
covariance tubes. Furthermore we visualize the early jumps after passage near a folded 
node in phase space for a model system with global returns. We also compute the proba- 
bility density of escaping trajectories on a cross-section for this example. SDKs have been 
integrated numerically by a standard Euler-Maruyama scheme [HS]. Deterministic 
solutions have been computed using a stiff ODE solver [HTJ 156] . 



7.1 Covariance Tubes 

The maximal canards and their associated covariance tubes can be computed. Fig ure [7| 
shows an example for these computations where we used the blown-up normal form (6.4) 
with /i = 0.08. 



First we compute the deterministic slow manifolds and by forward respectively 
backward integration (see [32]) up to the section 

S° = {{x,y,z) G : 2 = 0} . 

The primary and secondary maximal canards have been computed as intersections of the 
slow manifolds fl C^'; see also Section |3] and [531 EH]. The resulting maximal canards 
(thick green curves) are shown in Figure ITl We also computed a sample path (thin red 
curves) for each maximal canard starting at the same point as the maximal canard with 
Zq = —1; the noise values were fixed at a = 0.008 = a'. The tubes defined by the co- 
variance are shown in grey and have been computed using integration of the covariance 
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Figure 7: Computation of the canards and covariance tubes near a folded node in System 



(6.4). A detailed description of the figure can be found in Section 7.1 



differential equation (|6.8l); the section is drawn in yellow for better orientation. 
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Figure [7[a) shows a side view that illustrates how the different primary canards 7f 
and secondary canards 7^ are organized with respect to z. We have only started to draw 
the covariance tubes Bir) with = 0.02 a bit beyond the initial values at Xq = 1. It 
is clearly visible how the canards and their tubes are attracted towards the weak canard 
and then start to rotate around it. Figure [7](b) shows a front view towards the section S°. 
This view shows nicely how the tubes grow with increasing z-values and that the ellipses 
defined by the covariance matrix are indeed close to circular near z = 0. Furthermore we 
can see how the canards are organized on the maximal canard tubes for ^^'^'^ do not 
overlap while all other tubes overlap near the weak canard. Figure [7|^c) shows a zoom 
that illustrates the twisting and also shows how the sample paths are indeed "trapped" 
inside the covariance tubes with very high probability. 



7.2 Early Jumps 

To visualize the effect of early jumps, we consider a folded node with global returns given 
by 

dy = [-(^ + l)a; - z] ds + o' m?^ , (7-1) 
dz = [| + ax + 6x^] ds , 

which is a slight modification of a model system for folded-node MMOs [29] . The critical 
manifold is cubic-shaped (or S-shaped) and given by 

Co = {{x,y,z) eR^ -.y = x^ + x^} = Cq'^ U L_ U U L+ U C"^'+ 

where Cq'^ = Cq n {x < -2/3}, Q = Co n {-2/3 < x < 0}, C^'^ = Co n {x > 0}, 
L_ = Co n {x = —2/3} and L+ = Co fl {x = 0}. The parameters (a, 6) G help to 
adjust the global return mechanism. If a,b are 0{1) then they do not influence the local 
behaviour of a folded node at the origin (x,y,z) = (0,0,0). Figure [s] shows the effect of 
early jumps after passage through a folded-node region. Parameters for the simulation 
are: 

e = 0.01, = 0.143, a = 0.2, 6 = -1.1 , a = 0.005 , a' = . (7.2) 



In Figure [8[a) a deterministic trajectory (thick blue curve) has been computed for a = 
0. Then an SDK sample path for ( 7.1[ ) has been started on a point (green dot) of the 



deterministic solution and integrated forward. We define a cross-section 

:= {{x,y,z) G : X = -0.3} . 

Escapes of sample paths from the folded-node region are recorded on S"^. The next 
two returns are also shown as points (violet and red) on the cross-section S"'. It is clearly 
visible from Figure |8](a) how in this realization the SDK sample path jumps before the 
deterministic solution. Note that this causes the path to get re-injected into the folded- 
node region after a large excursion at a point slightly different from the deterministic 
solution. Hence the global return mechanism can potentially act as a control mechanism 
for the noise. To investigate the early jumps further we show in Figure [8|b) the probability 
density 

p{y, z) on E-^ . 
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Figure 8: (a) Phase space plot of simulations for (7.1). The deterministic solution (blue) 
and a stochastic sample path (black) are shown; parameters are given by (7.2) and for 
the deterministic solution we have a = = a'. Intersections with the cross-section E"' 
(yellow) are shown as thick dots; the SDE sample path is started on the deterministic 
solution and on S'^ (green dot). The next two intersections are shown as well (violet and 
red dots), (b) Probability density p{y,z) for 4000 sample paths on E"^ of sample path 
escapes from the folded node. The deterministic point-mass density is indicated as a bar 
Pdet (black). 



The density has been computed by recording the intersections with S after passage 
through the folded node for 4000 sample paths that have been integrated for a time 
s e [0, 20]. The corresponding deterministic point measure pdet has been indicated as well. 
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The density p{y, z) clearly shows that paths are likely to jump before the deterministic 
solution if we consider the z-coordinate distance from the folded node. We also see that the 
density p{y, z) is quite concentrated and shows a multi-modal structure. This structure 
can be explained from the fact that sample paths exit early but between different exit 
points they can make additional deterministic small oscillations. The possible different 
numbers of these oscillations correspond to the different maxima of p{y, z). 



8 Final Remarks 

In Section [6] we stated our results on the relation between the noise level, the parame- 
ter /i and the typical spreading of sample paths. Note that we proved and stated our 
results in blown-up (or re-scaled) coordinates removing the e-dependence. In particular, 
we worked in a neighbourhood of the folded node that is of size 0{y/e) in original co- 
ordinates. To obtain the results in original coordinates one has to apply a blow-down 
transformation. First, we replace {x,y, z,a,(T') by {x,y,z,a,a') (recall: we dropped the 
overbars for notational convenience), cf. (4.1) and ( 6.3[ ). Then the identity 



{x,y,z,a,a',i2) = {e ^/^x,e ^y,e ^/^z,e ^^V, e ^/V,/i) 
provides the required blow-down transformation. This implies, e.g., that the interaction 



of canards in Corollary 6.3 is given, in original coordinates, by relations of the form 



or that sample paths are likely to escape for 



z > 



log cr - ^ log e 



as shown in Theorem 6.4 Obviously one also has to translate the assumptions in a similar 
way, e.g., a <^ /i^^^ becomes cr <^ {^e)^^*. 

Another important point is that we focused on the detailed analysis near the folded 
node and did not consider different types of global return mechanisms. In Figure |9] we 
show the interaction between global returns and noise-induced early jumps on oscillatory 
patterns. The simulation in Figure |9] has been carried out using the model system (7.1) 
with parameter values 

e = 0.01, /i = 0.029, a = -0.1, b = -0.5 , cr = 0.005 , a' = 0.005 . (8.1) 

The deterministic solution in Figure |9] is an MMO with pattern 2^ where the number of 
small oscillations is difficult to count from the numerical results as /i is already extremely 
small. The SDE sample paths jump significantly earlier to C°'~ than the deterministic 
solution as expected from our results. However, we see that the deterministic solution 
makes an additional large-amplitude oscillation (LAO) given by the passage 

C"'" jump near L_ -> C"'+ jump near L+ C"'~ . 

In particular, the early jumps change the number of LAOs in the MMO pattern from 
L = 2 to L = 1. The noise also influences the number of SAOs but the crucial point is 
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Figure 9: Simulation for the model system (7.1) with parameter values (8.1). The upper 



plot shows a projection of a deterministic solution (black) and a stochastic sample path 
(red) into the {x, 2;)-plane. The lower two plots show the associated time series. 



that it can also have a global effect. Hence we end up with a stochastic- resonance-type 
mechanism for MMO patterns. 



Furthermore, one could think about extending our results to capture the effect of a 
global return map on the escape density after the passage through a folded-node region; 
see also the brief discussion in Section l7.2l Let S"^ denote a cross-section on which we 
record the escape from the folded-node region and let denote a cross-section to the 
deterministic flow slightly before the re-entry to the folded-node region; see also [79] . 
Then we have a global return map induced by the deterministic flow 

M : ^ . (8.2) 

If we are given a probability density p on S'^, we can then consider the induced density 
p o M~^/|detDM o M~^\ on T^^. This should allow us to calculate a distribution for 
different MMO patterns, i.e., we can hope to assign a probability to each combination 
after one return from E"^ to S"^. Although this approach seems possible it is beyond the 
local analysis we focused on here. A solution of this problem crucially depends on the 
form of the global returns which are described by the map M. 
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A Proof of Theorem 4.3 (Canonical Form) 



We are going to need the following lemma for the proof of Theorem |4.3[ 

Lemma A.l. Consider two non- autonomous vector fields F,G : M x — )• M.^ for 
{z,X) G M X M^. Suppose that on an open set T>, containing (0,Xo), both are continuous 
in z and in X. Furthermore, suppose for all {z,X) in D we have 

\\F{z,X) -G{z,X)\\ < 12 . 



Let K be a Lipschitz constant for F in X, which is uniform in z. Suppose X{z) and Y{z) 
solve ^ = Fiz.X) and ^ = G{z,Y) with X(0) = Xq = Y{0), then 



\X{z)-Yiz)\\^^{e^\^\~l) . 



(A.l) 



Proof. A direct Gronwall-lemma argument suffices; for details see [60], p. 399-400. □ 



The proof of Theorem 4.3 proceeds in several steps that aim at bringing Equation (4.8 ) 



into diagonal form [T71 126] . 
Proof of Theorem As a first step we apply a rescaling 



u = exp 



2/i 



Tr(A(s)) ds 



where u^^^ = (w^^^Ug^^) are new coordinates. This yields 



II— — = Ao{z)u^'^', with Ao{z) 



2z 2 
-2(/i+l) -2z 



(A.2) 



where now Ti{Aq{z)) = 0. Therefore the principal solution has determinant 1 and is area 



preserving. We are going to show that the solution of (A.2) is a rotation up to a small 



error using a sequence of z-dependent coordinate transformations Sj{z). We set 

-z—iu]{z) —z+iuj{z)^ 



= So{z)u^^\ where So{z) 



i+fj. 
1 



1 



Then Sq ^AqSq is diagonal with entries ±2ia;(z). We get 

dM(i) 



dz 



Ai{z)u 



(1) 



where the new matrix Ai is given by 



Ai{z) = Sq ^AqSq - iiSq 



i-uj'{z) 



dz 



,, -iW(^) 
2uj{z) 



,i_±Az) 

2aj(z) 

—2iuj{z) + /i- 



'i-Uj'(z) 



2uj{z) 



As above we want to have a zero trace so we compute 



2^ 



Tr(Ai(s))ds 



2 Jo uj{s) 
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ds 



\ogu{z) + const . 



Hence we consider the scaling 



u 



= ^i(^) 



u 



(2) 



u 



(2) 



which yields the equation 

d«(2) 



with Aoiz] 



where overbar denotes complex conjugate and 

/i 



0J2{z) = 2uj{z) + 



2u{z) 



i+uj'{z) 



2u{z) 



(A.3) 



The next transformation 
yields 



u 



u 



(3) 



IX— = A^{z)u'^^\ 



A:,(z) 



S2' 



A2S2 - /i 



d^ 

dz 



Instead of a given transformation we now impose the form of the matrices 

1 jjiv{z) 



S.iz) 



flV[Z) 



Asiz) 



Piiz) 
Piiz] 



(A.4) 



(A.5) 



so that ^2 = ld+0{p) and ^3 is diagonal. Substituting (A.5) into (A.4) leads to the 
equations 

= iu2{z) + fi^p2 {z)v - pi {z) 

P^z = P2{z) -iUJ2{z)v -Vpi{z) 

and their complex conjugates. The first equation determines pi{z), and thus the second 
one becomes 

-2iu2{z)v - p^P2{zy + P2{z) . (A.6) 



dz 



If we can show that (A.6) has a bounded solution for times \z\ < 1 then our prescribed 
coordinate change 5*2 exists. Now let a{z) = 2z + Re pi{z) = 2z + 0{p^) and define Zq) 

i r 



by (4.24). Then a last transformation 

«(3) = S^{Z)U 



-«(^,o)/m 



1 



1 + iVl i 



u 



(A.7) 



will bring the equation into canonical form (4.21), with w{z) = Impi(2;) = 2uj{z) + 0{p). 
Composing all transformations, we get u = S{z)u, where S{z) = e^ 6*0 (2) Si (z) 5*2 (z) 6*3 (2;) 
is indeed of the form (4.20). 



To prove the existence of bounded solutions for (A.6) note that P2iz) and cu2{z) are 
bounded away from zero and that they are also have bounded norms for \z\ < 1. Now set 

N 21^2(2) , P2iz) 

F{z,v) := V -\ , G{z,v) :-- 



2iu}2{z) - r ^ 2 , P2{z) 

-V - PP2[Z)V + 



p p p p 

Note that F is a vector field with bounded solutions on an (9(l)-time scale. Considering 
-F, G as real planar vector fields on M^, we can apply Lemma A.l to conclude that (A.6) 
admits bounded solutions on time intervals of length 1. □ 
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B Proof of Theorem 6.1 (Covariance Matrix) 



In this Appendix, we discuss the equation 



?0\T 



dz 



(B.l) 



describing the evolution of V{z) = {Vij{z)}i j ^^12} := o"^^ Cov(2;), the covariance matrix of 
the hnearized variational equation around a deterministic solution of (6.4). Before proving 
Theorem 6.1 on the small-// asymptotics of the solutions, we provide different approaches 
yielding information on the behaviour of V{z). A formal method based on iterative 
computations of a slow manifold is developed in Section B.l to understand the asymptotics 
of (B.l ). In Section B.2 we provide rigorous bounds using a Lyapunov function. In Section 



B.3 we refine the previous results by a transformation to real canonical form and results 



about delayed Hopf bifurcation, thereby proving the main theorem. 
Proposition B.l. Let 

vi{z) = Vuiz) 



v{z) 




where 



V2{z) = V22{z) 

y^[,z) = Vi2{z) = 1/21(2;) 

denote the two variances and the covariance. Then v{z) satisfies the ODE 

2 



(B.2) 



dv 

d^ 





v + 



■■.B{z) 




(B.3) 



where we have abbreviated x'^^^{z) =: x{z). 



Proof. Using the definitions from (6.6) and equation (B.l) we get 
dV 



dz 



A{xf')V + VA{xfY + (F°)(F°)^ 

_4^det _2(^ 



v + v 



_4^det 2 

-2(/x + l) 

AV12 - SVuxf^ + 2 
. 21^22 -4\/i2a;f*-2\/n(l + /i) 

Therefore the result follows. 



r 







V2 

V2p 

2F22-4\/i2xf*-2yn(l + /i) 
-4\/i2(l+/i) + 2p2 



V2 
V2p 



Our goal is to analyze (B.3) for a given maximal canard solution x{z) 
B{z) has eigenvalues 



□ 



Observe that 



-Ax{z), -Ax{z) ±Ai^Jl- x{zy 
We assume that Zq, Zi are chosen so that 

-l<2;o<0<2;i<l and 1 - x(2;)^ + /i > 



4x(z) ±4:iujo{-x{z)) . 



Mz e [zq, zi] 



(B.4) 



In particular, the assumptions (B.4) are satisfied for any maximal canard solution ap- 
proaching the folded-node region from the slow manifold C° for some 2:0 < < 2;i of order 
1 and /i sufficiently small. 
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B.l Iteration and Asymptotics 

Notation B.2. Henceforth, we write x{z,fi) x y[z,fi) if 

c-y{z,n) ^ x{z,n) ^ c+yiyz.jj) 

holds for all z with positive constants c± independent of z and /x. 



(B.5) 



A formal derivation for the asymptotics as /i — > for (B.3) can be carried out using 



an iterative scheme [901 [IT] . We set v{z) = Vq{z) + Vi{z) where 

1 



V*{z) := -B{z)-'E 



A{l + iJ,)x{z) 



1 + /i + 
1 + + {1+Ax{zf + /i)p2 
2p'^x{z) 



(B.6) 



defines the critical manifold for (B.3) when viewed as a slowly time-dependent system. 
We get 

^^^ = Biz)V^ + ^xE^iz), E,{z) 



dz 



The same change procedure also works for any n ^ 1 by setting 

n 
3=0 

Then Vn+i{z) satisfies the equation 

dK+i 



(B.7) 



dz 



(B.8) 



where V* and En are given inductively by 

v:^,{z) = B{z)-'j-v:{z\ 



Er 



dz 



[B{z)-'En{z)] 



Remark: Observe that (B.7) is the asymptotic expansion for the slow manifold of (B.3). 
The iterative scheme we use here is very convenient for slowly time-dependent systems. 
Many other methods to calculate slow manifolds for general fast-slow systems have been 
explored; see |110[ 1111] and references therein. 

Proposition B.3. Assume that the deterministic maximal canard solution x{z) satisfies 
(B.4). Then the asymptotic expansion (B.7) of v{z) forn ^ has components of order 



^-v:jz)^fxw:Jz) 



p" 



|2n+l ' 



|2n 



(B.9) 



Proof. First observe that the symmetry (3.12) implies that x(0) = for any maximal 
canard. Using this fact and the form of the slow flow (3.10), we find that x{z) must have 
a Taylor expansion with non- vanishing linear term, i.e.. 



x{z) = XiZ + X2z'^ + 



(B.IO) 
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with Xi < 0. The proof of (B.9) then proceeds by induction as follows: The base step 
n = holds by formulas (B.6) and (B.IO). The induction step from n — 1 to n follows 
from direct differentiation 

2n- 1 

^2n ' 



dz \ z 



y2n-l 



and the calculation of B{z) ^ 



Biz) 



8(1+At)x{2) 
1+M 



+ 



\ 



8x(2) %x{z) ' 2(1+Af) 





1 



4(1+m) 



0\ 

1 

2 



almost immediately; we just have to observe the block structure of Bi^z 



,-1 



□ 



Proposition B.3 is a formal asymptotic result. The asymptotic series (B.7) becomes 



"disordered" for \z\ = 0{y/Ji), because in this case all the terms for the coordinates vi 
and V2 are of order l/y^, while all terms for ^3 are of order 1. Therefore we conjecture 
that 

„, = 0(-l=), ". = 0(-l=), „3 = 0(1), (B.11) 

for —y/JI < Zo ^ z ^ y/JI. 

B.2 Lyapunov Function 

The results in this section are not as sharp as the results obtained by coordinate changes 



in Section |B.3[ but they are obtained by a completely different technique which is of 
interest on its own in the context of folded nodes. Next we are going to establish an 



auxiliary result needed below in the proof of Proposition B.5 



Lemma B.4. Consider the linear non- autonomous differential equation on M given by 

(B.12) 



ijL— — = kizX + k2 



dz 



(-^) 



2n 



where ki^2 = (if^G two positive constants, ^ > 0, and either n = and z ^ zq or 

n ^ 1 and zq ^ z < 0. Then 

|-(2n+l) for Zo + O{fl\\ogfl\) ^ z ^ 

X(z) X <! for ^z<Oifn^l, 

for — JJi ^ -2 ^ JJi if n = . 



-1/2 



Proof. The solution of (B.12) can be written as 



where 



X{z) = X{zo) e'^^(^'-^o)/2'^ +k2ln{z 

1 



(B.13) 



/i 



n-1 



^fcl(22-i2)/2M 



20 



-t) 



2n 



dt . 
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For z e [zo,— a/zu] the leading-order asymptotics of (B.13) is given by /„(2). Using 
integration by parts, we get an upper bound 



fx" 



n-1 



1 



(-t)2"+l 



kit 



' d 

dt 



l2n- 



(B.14) 



The lower bound follows by inserting the upper bound for In+i{z) in (B.14). Note that 
here the condition z ^ zo + (9(/i|log is needed to make the term e'^i(^^~^o)/2M small. Let 
us also remark that, in particular, we showed X{—y/JI) x 1/ ^JJi. 

Finally, to describe the behaviour for —JJi < ,2 < 0, we replace Zq by —JJl in (B.13). 



Then all exponential terms are of order 1, and the integral can be estimated directly. □ 



Proposition B.5. Suppose (B.4) holds, and let x{z) he a maximal canard solution. Then 
solutions to the variational equation (B.3) remain hounded hy 0{l/{\z\ + ^/Jl)) for zq ^ 

z ^ y/JI. 



Proof. Throughout the proof we are going to introduce several positive constants cj = 
0{1) for j G N, whose actual values do not influence the asymptotic result. As a first 
step, we want to find a symmetric matrix M{z) such that 



B{z fM{z) + M{z)B{z) = -x{z) Id 



(B.15) 



This can simply be accomplished by solving the six algebraic equations (B.15). We find 
that 



M{z)\ 



1 



ti=0 



64(1 + 3x2 



^7 + 12x2 1 + 12x2 
1 + 12x2 7 + 64x2 + 48x^ 
-12x -16x(l + 3x2) 



-12x 
-16x(l + 3x2^ 
4(3 + 18x2^ 



where we have abbreviated x 



xiz) 



Since (B.IO) holds for maximal canards, it is 



straightforward to check that the matrix M{z) is positive definite, uniformly in n and z. 
Therefore it defines a family of quadratic forms 



r„ := V^M{z)Vn 



(B.16) 



where Vn is defined by (B.7). The quadratic form (B.16) satisfies 

for some constants Ci,C2 > 0. Essentially Yn will act as a Lyapunov function to bound 
To show this, we compute the derivative. Using (B.15) and (B.8), we get 



+ t^V^ 



AM. 



■Vn 



AM 



x{z)VjVn + l^Vn^Vn + /i" [En{zf M{z)Vn + M{z)En{z)\ . (B.17) 
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Since \\Er, 



and II 'i^'^ I 



is bounded, we can find constants 03,04 > such that 



(B.17) imphes 



dK 



dz 



Setting F„ = Z^, we find that the last inequahty is equivalent to 



2n 



(B.18) 



Using (B.IO) and Lemma B.4 we obtain that for Zq of order — 1, 



ZJz) ^ ct 



for zq ^ z < 



(B.19) 



Since Zn is equivalent to this shows that (B.8) is indeed an asymptotic expansion 

in powers of fi/z^ for z ^ —y/Jj-, and in particular all components of v{—y/fl) are of order 
To complete the proof up to time z = ^/Jl, we simply apply (B.19 ) in the particular 
case n = (that is, for Vq = f and Eq = E). 



□ 



tion 



In view of our conjecture (B.ll), the bound on the covariance provided by Proposi- 



B.5 is not sharp since we have not yet shown that = 0{1). 



B.3 Delayed Hopf Bifurcation 

To obtain a sharp bound on the covariance we consider a similar coordinate change idea 



as in Section 4^ This procedure will give a variational equation for the covariance that 
has desirable symmetry properties. 

Lemma B.6. There exists a linear coordinate change = S{z)(^ transforming the lin- 



earized SDE (6.6) into 



dCl = -A{z)~Cdz 



Fiz) m... 



(B.20) 



where A{z^ is in canonical form 

A{z) 



a(^z) w{z) 
-uj(z) a(z) 



(B.21] 



with a{z) = —2x{z) + 0{fi^) and w{z) = 2uj{z) + 0{fi). The matrix F{z) is positive 
definite, with eigenvalues bounded below and above uniformly in z. 

Proof. It suffices to apply the coordinate change C° = S{z)C^ constructed in the proof 



of Theorem 4.3| (with an obvious modification due to the fact that x{z) is not necessar- 



ily given by the weak canard). The new diffusion coefficient is then given by F{z) = 
5(z)-iF0. □ 



A computation analogous to the one in the proof of Proposition B.l then yields 
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Lemma B.7. The covariance matrix of is given by a'^V{z), where the matrix elements 
of V{z) satisfy the system 



2a{z) 2tu{z) 

2a{z) -2w{z) \v + E 
—w{z) zu{z) 2a(z) 



(B.22) 



--■.B{z) 



where E vector in with 0{1)- components. 



It is already apparent from the form of (B.22) that the analysis of the variational 



equation simplifies. We can now prove Theorem |6.1[ which we restate as follows for 
convenience. 



Theorem B.8 (Theorem 6.1). Suppose (B.4) holds. Then the solution v = (^1,^2,^3) 

-> 

^3 = 0{1), 



for (B.3) satisfies the following asymptotics as /i — > 
1 1 



V2 



\z\ + y/i \z\ + y/JI 

for zo + C(/i|log/i|) 



{v, - V2) = 0{1) 



(B.23) 



equations of (B.22) we get 



Proof. We work in the coordinates provided by Lemma |B.7[ Summing the first two 

(B.24) 



/i— (wi + V2) = 2a{z){vi + V2) + 61 + 62 
az 



and we already know from Proposition B.5| (resp. Lemma B.4) that this yields 

{.Vi+V2){z) 



\z\ + ^ 



for 2:0 + (9 (/i I log /i I) ^ -2 ^ ^/Ji. The difference of the first two equations in (B.22) and the 
third equation can be combined as 



/i 



dz 



Vi - V2 

V3 



2a{z) Aw{z] 
-wiz) 2a(z) 



V3 



ei - 62 
63 



(B.25) 



Considering (B.25) as a fast-slow system with slow variable z we find that the critical 



manifold is given by the equation 



Vi -V2) 



2a{z) Aw{z] 
-w(z) 2a(z) 



-1 



ei - 62 
63 



which is of order 1. Observe that (B.25) undergoes a delayed (or dynamic) Hopf bifurca- 
tion at z = 0. Thus Neishtadt's theorem on delayed Hopf bifurcations |90] applies, and 
shows that solutions of the variational equation satisfy 

(V, - V2)iz) = {Vi - V2nz) + 0(/i), Vsiz) = V;{Z) + 0(/i) 



for 2o+C(/i|log/i|) ^ z ^ C(l). Now the result (B.23) follows from ^(2) = S{z)V{z)S{z 



by writing V{z) as the sum of a leading term proportional to the identity matrix and a 
remainder of order 1. □ 
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C Proof of Theorem 6.2 (Staying in covariance tubes) 



Applying the transformation of Lemma |B.6| to the nonhnear equation (6.12) and dropping 
the tildes yields the system 



dC. = - [A{zK, + 6(C., z)] dz + -^F{z) dW, 



where 



A{z) 



a{z) w{z) 
-voiz) aiz) 



(C.l) 



(C.2) 



and h{C,,z) = 0(||C|| )• The solution of (C.l) with initial condition ^2^, = can be written 
as 

Cz = ^ f Uiz, s)F{s) dW, + - r U{z, s)b{C, s) ds =: C° + Cl , (C.3) 

where U{z, s) denotes the principal solution of the time-dependent linear system fi( = 
A{z)(. Owing to the particular form of A{z), we have the explicit expression 



jjt^ g\ =Q-c.{z,s)/,, f cos{<^{z,s)/fi) sm{<^{z,s)/fi)\ 
\-sm{Lp{z,s)/fi) cos(v?(z, s)//i)y 



where 



a[z,s] 



-a(u) du 



ip{z, s) 



w(u) du 



(C.4) 



(C.5) 



Note in particular that since —a{z) x x{z) x —z near z = 0, we have a{z, s) x — z"^. 

For a two-by-two matrix M, let ||M|| denote its L^-operator norm, i.e., ||Mf is the 
largest eigenvalue of MM'^ . 



Lemma C.l. Let 



Then 



Q{z) = \ f\\U{z,s)\\ds 

H- Jzo 



eiz) = o 



\z\ + ^ 
for all Zq ^ z ^ ^JJi. 

Proof. Since [/(z, s)'^ = Id, we have 



Q[z) = - / e-"("'^)/'^ ds 

Jzo 



(C.6) 
(C.7) 



and the result follows from Lemma [B. 41 

The next lemma provides bounds on the norms of V{z) and V{z)'^. 
Lemma C.2. Let 

K^{zr=\\V{z)\\ , K_{zr = \\V{z)-'\\ . 

Then 



□ 



K+{z) 
for all zq ^ z ^ y/JI. 



O 



1 



1^1 + 



K^{zf = 0{\z\ + ^) 



(C.9) 
(C.IO) 
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Proof. Note that 



viz)viz: 



Vi V3 
V3 V2 



{Vi + V2)V3 v\ + vl 



has eigenvalues given by 



(C.ll) 



(C.12) 



The larger eigenvalue is equal to -^'+(2;)^, while the smaller one is equal to K^i^z) ^. The 
result thus follows from the bounds obtained in Theorem IB.8[ □ 



We can now prove a local version of Theorem 6.2, on a small interval C [zq? A- 



Proposition C.3. Fix times zo ^ s < t ^ z such that a{t,s) ^ C>{fi). Then, for all fi 
and a small enough, for any < 7 < 1, 



P<; sup {Cu,V{u)-\u)>r 

,2 



1 



1-7 



exp 



7r 
'2^ 



l-0{i\s\ + ^ 



t-s 

I 

/i 



. (C.13) 



Proof. The proof is adapted from [211 Section 5.1.2], and we use almost the same notations 
as there. Let T„ = U{s,u)(u- Then 



{Cu,V{ur'Q = {T^,U{u,sfV{ur'Uiu,s)T^) = \\Qs{u)TA 



(C.14) 



We can decompose T„ = T° + T^, where 



T° = ^ / U{s,v)F{v)dW,, 
1 

Tl=- U{s,v)b{C,v)dv. 



(C.15) 



The process is a Gaussian martingale. Lemma 5.1.8 in ^24j can thus be applied and 
provides the bound 



pjsup ||Q,(t)T?|| ^i?oj ^ -^exp(- 



^2a2 



(C.16) 



For this bound to be useful, we need to show that Qs{u) and Qs{t) are close to each other. 
Observe that 

^i^Qsiuy = f/(s, u)F{u)F{ufUis, uf , (C.17) 
au 

as a consequence of the definition of U{s,u) and the differential equation satisfied by 
V{u). Integrating from u to t and multiplying on the left by Qsiu)"^, we get 



1 /■* 

QsiuyQsit)-^-ld = Qsiuy- / 



U{s,v)F{v)F{vyU{s,vy dv . 



(CIS) 
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Now ||f/(s,i;)|| ^ 0{1) owing to the assumption a{t,s) = 0{fi). Thus the integral has 
order t — u ^ t — s. Furthermore, 



\\Qs{u)'\\ ^ \\U{u,s)fK^iuf = 0{K^{uf) . 
As a consequence, we get 

d — s 



Id+O K(s 



(C.19) 



(C.20) 



Thus there exists an i? = r[l — 0{K_{sY{t — s)/ jj)] such that 

sup (C„,F(n)-iQ sup ||Q,(t)T„|| ^ i?| . (C.21) 

For any decomposition R = Rq + Ri with Rq, Ri > 0, we can bound the above probabihty 
by Pq + Pi, where 



Pi = P<i sup \\Qs{t)Tl\\^ Ri 



i = 0,1 . 



(C.22) 



Po has already been estimated in (C.16). Now we want to choose Ri in such a way that 



Pi = 0. For any u ^ t A tq(j.)^ we have 

1 r 

\\Qs{t)^l\\ ^ const sup K_{u)- / \\U{s,v)h{Q\\dv 



^ const sup K^{u)Q(u) sup 

^ const sup K^{u)Q{u) sup fr+(f)^r^ 

u£[s,t] v£[zo,t] 

^ const {y/Jl + \t\y^^^ . (C.23) 

We can thus achieve Pi = by simply choosing Pi as a sufficiently large constant times 
{y/JI + |t|)^^/^r^. This determines Rq, and the result then follows from (C.16). □ 



We can now complete the proof of Theorem 6.2 which will follow directly from 



Theorem C.4. There exist constants Ao,ro,/io > such that for all < A < Aq, all 

< a < r < r^jj?/^ , < jj < and all < 'j < 1, 

P{rew <z}^ C'+(z,;2o)exp|-7^[l - 0(A) - ^(r/x-^/^)] | (C.24) 



holds for all z ^ y/Jl, where 



C+{z,zo) 



const 
(l-7)A/i 



* ds 



(C.25) 
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Proof. Let Zq = Sq < Si < ■ ■ ■ < = z he a. partition of [2:0, z]. Then 



N 



P{rB(.)<s}^^Pfc 



k=l 



where 



Pfc=p( sup {Cu,Viu)-\^)^A 



(C.26) 



(C.27) 



can be estimated by Proposition |C.3 We want to choose the partition in such a way that 
the error terms in Pk are bounded uniformly in k. A convenient choice is to define the Sk 
by 

a(sfc+i, Sk) = Ajj, when Sk+i < , 
Sfc+i - Sfc = when ^ y/i 



(C.28) 



Using the fact that s) x + — s) and applying Proposition C.3, one indeed checks 
that 

„2 



,...,iV, 



(C.29) 



where the error terms are uniform in k. It remains to estimate the number of elements 
of the partition, which will give the prefactor C+. In the case z ^ — a//^, we simply have 

Ol{z,Zq) 



a{z, zq) = NAn 



N 



In the case —^/Jl ^ z ^ ^/Jl, we have 

' a{z,zo) ~\ p- (-yZ/i) 

and the result follows from the fact that a{z, —^/Ji) x ^/Ji{z — {—^JJi)). 



(C.30) 

(C.31) 

□ 



Theorem 6^ is just a reformulation of this result, in which we have chosen 7 
1 - a^/r\ 



D Proof of Theorem 6.4 (Early Jumps) 



We consider again the equation for the difference C,z between stochastic sample paths and 
a deterministic reference solution, this time given by the weak canard. In canonical form, 
we have 

dC. = - [A{z)Cz + 6(C., z)] dz + -^F{z) diy, , (D.l) 

where 



and b{(,z) = 0(||C|| )• The proof of Theorem 6.4 is split into several parts. In Sub- 



section D.l , we show that sample paths are likely to leave a neighbourhood of order 
slightly (that is, logarithmi cally) larger than c/y/z of the weak canard in a time z of 



order a/ /i|log/x|. Subsection D.2 analyses the dynamics in a larger neighbourhood of the 



weak canard, in which the drift term dominates. Subsection D.3 combines the two results 
to prove the main theorem. 
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D.l Diffusion-Dominated Escape 

We assume from now on that a <^ /^^''^, because otherwise stochastic sample paths are 
no longer localised near deterministic solutions when z = ^Ji. We define the set 

S{h) = {{C,,zYz^^,\\a<hp{z)} , (D.3) 

where 

The following result is an adaptation of pT| Proposition 4.7] to the two-dimensional case. 
Proposition D.l. Let h,!/ > 0, with v of order 1, satisfy the conditions 



a 
h 



^ Co and 



2+u 



log 1 + Z/ + 



a' 



1/2 



3/4 



a 



(D.5) 



/or some Co,Ci > 0. //cq and Ci are small enough, then there exist T > and C(z/) > 
such that for any {zq, Co) G S{h) with zq < T, 



h 



a 



2u 



exp< —K^u) 



(D.6) 



holds for all ^JJi ^ zq ^ z ^ T , where 



l-O 



v\og{h/a) 



(D.7) 



We shall choose the value of the parameter u later on, while h will be taken of the 



form h = CO" I logo" I . Then Condition (D.5) reduces to 

a|log(Tp+Vlog|log^|^0(/^'/') , 



(D.8) 



which is slightly stronger than requiring a <^ ji^^^. The exponent k(z/) in (D.7) becomes 
optimal in the limit v — )■ oo, but Condition (D.8) becomes more stringent as v grows 
large. Note that we have to choose a finite v anyhow. 

Proof of Proposition D.l . Let a{t, s) = a{u) du. We define a partition zq = Sq < si < 
■■■ < sn = z oi [zq, z] by 



= A/i , /c = 1, . . . ,iV - 1 , 
where A > will be chosen later. The Markov property implies that 



(D.9) 



where 



\ = sup P^'^ 

C: IICII^M^fc-l) 



k=l 



IICsll ^ , 

sup — — < n 



(D.IO) 



(D.ll) 
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We shall derive a uniform bound q{A) for all Pk, 1 ^ k ^ N — 1. Then (D.IO) and the 

a{z,zo) logg(A)-i 



definition (D.9) of the partition imply 



/i 



(D.12) 



and the result will follow from an appropriate choice of A. 

The process C,s starting at time s^-i in C, can be decomposed as = Cs'" + with 



U{s,u)F{u) dWu 



c, 



k,l 



fx 



U{s,u)b{Cu,u) du 



(D.13) 



Sk-l 



where the principal solution U{s,u) of the linear system is given in (C.4). For any de- 
composition h = Hq — Hi, we have 



Pk ^ sup [Pk,o{C, Ho) + Pfc,i(C, Hi)] 



(D.14) 



where 



Pk,oiC,Ho)=f^'''''{ sup 



fc,0| 



<Hn 



> 



Pk,iiC,Hi)=F^'' 



fc,ll 



sup . 



lie 



^ Hi, sup 



< h 



We start by bounding Pk o{(, Hq), using the end-point estimate 

nH^p{sk)^ 



(27r)2detCov(0?) 



(D.15) 



(D.16) 



The last inequality follows from the fact that the random variable '° is Gaussian, and 



Sk 



we have bounded its density by the normalizing constant. We denote the diagonal matrix 
elements of = cr^^ Cov{Cgf) by Vi and V2, and the off-diagonal element by ^3. Then 



detCoyi& = a\viV2-vl)=a' 



1 



i{TTV)'-ivi-V2)')-vl 



(D.17) 



As already remarked in the proof of Theorem B.8 , the quantities vi — f 2 and f 3 remain of 
order 1 up to some 2; = T of order 1, as a consequence of Neishtadt's result on delayed Hopf 
bifurcations. In order to estimate the trace Tr V, we use the fact that p{z) is decreasing 
in (0,T] for T small enough, owing to the fact that Tr {F{z)F{z)'^) is bounded below 
by a positive constant, and has a derivative bounded in absolute value. Thus we have, 
cf. (|B^, 



TtV 



/i 

„2A 



^2a{sk,u)/fj. 



Tr(F(M)F(M)^) du 



Sk-l 

Au 

k-i 

2 r„2A 



^p(.,)^[e^^-l] . 



(D.18) 



44 



Substituting (D.18) in (D.17) and then in (D.16) yields 



l + O 



(j2 e2A_l 

Next we estimate Pk,i- We first obtain the bound 



(D.19) 



Aft, 



\\Uis,u)\\ \\biCu,u)\\du 



Sk-1 



< const- 



Sk-l 



< const P^''^'^' -^^ 



2s 



k-l 



(D.20) 



where we have used the fact that the function u i— )■ p{uY/2u is decreasing, and bounded 
the integral of {2u/ jj) e"'^^''^^^^ by e^. Using a Taylor expansion of and the definitions 
of p and of the partition, one finds 



^ 1 + 



S - Sk-1 



(D.21) 



for all s G [sk-i,Sk]. Together with (D.20), this implies 



Pis) 



< const 



p{sk~i) /— A . . /i vAe 



-VA e'^ ^ const 



Sk-1 



P 



3/4 



2 



which yields -Pfc,i(C) Hi) = 0- Substituting Hq = h + Hi in (D.19) thus yields 



Pk < g(A) := 



2 o2A _l 



e 



l + O 



1 



/i^Ae^ 



P 



3/4 



for A; = 1, . . . , — 1. Finally, we make the choice 



l + iy 



log 1 + Z/ + 



Bounding above amounts to bounding q{A) ^ below. For this we write 



g(A) 



-1 



a 



,2A 



1 -0(e 



-2A^ 0(e-4A) _ 

2\ 1+v 



/iv^e^ 



/i 



3/4 



1-C 



O 



P' 



h fh 
a 



3/4 



log 1 + Z/ + 



(D.22) 

(D.23) 
(D.24) 

(D.25) 

/^2xl/2 



Note that the error term 0{e is negligible and no longer appears in the last line. Now 
by Assumption (D.5), for cq and Ci small enough, we get 



a' 



(D.26) 



and 



logg(A)-i ^ 2z/ 



2 log 2 



A 



1 + (l + z/)log(l + z/ + /iVc^^) 



=: k{v) . 



(D.27) 

Note that k(z/) is indeed of the form (D.7). The result thus follows from (D.12). □ 
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D.2 Averaging and Drift-Dominated Escape 



We consider again Equation (D.l), but this time for slightly larger values of We start 
by transforming the system to polar coordinates. 

Lemma D.2. Consider a system of the form 



(D.28) 



•where Wz denotes a two-dimensional Wiener process, and and are row vectors of 
dimension 2. Then in polar coordinates = rcos(p,ri = rsimp), the system becomes 

dr = -fr{r, ip, z) dz + -^Fr{ip, z) dW^ , 
/i ^ 

XX (T X 

d'P = --f^{r,if,z)dz + ^-F^{^,z)dWz , (D.29) 



where the new and old diffusion coefficients are related via 

Fr = F^ COS ip + Fr^simp , 
F^ = F^ cos ip — F^smip , 

while the drift coefficients are given by 



(D.30) 



fr = COS !p + frjSm!p+ ^F^F^ , 



fif = fr, COS if- f^ sin if 



a 



-FrF^ . 



(D.3i; 



Proof. The formulas can be checked directly by applying Ito's formula to (D.29). 



□ 



Applying this result to (D.l), we obtain a system of the form 

.2- 



dr 



dip 



1 
/i 
1 



a{z)r + r'^briip, z) + O 
-w{z) + rb^{ip, z) + O 



a 



dz H Fr((p,z) dWz 



2\-\ 



dz+ ~F^{^,z) dW, 



(D.32) 



Note that the functions br and 6^ do not depend on r, owing to the fact that the nonlinear- 
ity in the original equation is homogeneous of degree 2. Another important observation 
is that the average of br (and b^) over ip is zero. This follows again from homogene- 
ity, combined with (D.31). This observation suggests to simplify (D.32) by an averaging 
transformation. 

Proposition D.3. There exists a function w{ip,z), which is bounded, smooth, and 2it- 
periodic in ip, such that f = r + r'^w{ip, z) satisfies the SDE 



dr = - \aiz)r + /3(r, ip, z)] dz + -^Fr(r, ip, z) dW^ , 



(D.33) 
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where 



Frif,cp,z)=Friip,z) + Oif) . (D.34) 
Proof. Using Ito's formula and the fact that r = f — f^w{ip, z) + (9(f^), we obtain 



dr 



+ 



dw dw 

a{z)f + ( 2zw{ip, z) + hr{(p, z) - w{z) — + /i— ) + 0{f^) + O 



a 



d(p dz 

(dw \ 
2w{^, z)Fr{^, z) + -Q^F^iv, z) \ + 0(f2) 



d^; 
(D.35) 



It is thus sufficient to show that the equation 



dw a{z) br{ip,z) 

dip w[z) zu{z) 



(D.36) 



solution of (D.36) can be written 



admits a bounded, 27r-periodic solution. Letting c = c{z) = a{z)/w{z), the general 

(D.37) 



w{(p,z) = w{0,z)e'''^ + 



Thus choosing 



,27rc /•27r 



w{0,z) 



1 - e 



2ttc 



-ap 



■Dj{z) 
br{ip,Z 



W[Z] 



dip 



the resulting w{ip, z) is indeed 27r-periodic as a function of (p. Finally note that 



lim ti?(0, z) = — 



27r 



U7(0) 



dip 



showing that w{ip, z) is also bounded as z — )■ 0. 
We now define the set 

V{ri) = {{f,ip,z): z^ y/Jl,f < r|^/z} . 
Then the nonlinear term P satisfies, on the set TZ = T>{ri) \ S{h), 



\f3{f,^,z)\ 



rz 



^ M - + 



—2 — 2 

r jir a 



z z r^z 



3/4 



|loga| 



(D.38) 

(D.39) 

□ 

(D.40) 
(D.41) 



for some constants M, M' . We can thus find, for any < k < 2, an 77 = ?7(k) such that 
(for sufficiently small a and //) 



a{z)r + /3(r, ip^ z) ^ Kzr 



holds in TZ. Thus we have 



dr = — \Kzr + /3(r, y9, z)l dz H -^r(^, V^, z) dWz 



(D.42) 



(D.43) 
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where (3{r, {p,z) ^ in 71, implying 



+ ^ r e-(^'-^'y^f^Fr{f,,ip,,s)dWs {BAA) 



holds as long as the process stays in TZ. 

The following proposition shows that the process is unlikely to stay in TZ for times 
significantly larger than (/i|log(j|)"^/^. 

Proposition D.4. There exists a constant ^2 > such that for any initial condition 

(fo, v^o, zo) G n, 



{t-ji > z} ^ 2expl -K2 



' fi\\oga\ 

Proof. We introduce a partition zq < zi < ■ ■ ■ < = z of [zq, z], given by 



^fc+i — zl= 7/i|loga| for ^ A; < = 
The Markov property implies that 

N-l 

p(fo,c,o,.o)|^^>^| ^Y[Pk 

where 



7/i|logcr| 



k=l 



Pk = sup Pkir,ip) , 

f,ip: (f,(p,Zfc)e7?. 



Inequality (D.44) (with zq replaced by z^) shows that 



a 



fz, + —M:, 



fz > 

holds for Zk ^ z ^ T-ji, where is the martingale 

= I e~^^''-''^y^^Fr{fs,iPs,s) dWs . 

J Zh 



It follows that 



(D.45) 



(D.46) 



(D.47) 



(D.48) 



(D.49) 



(D.50) 



(D.51) 



where we have used f ^ hp{zk) = ca\\oga\p{zk) and e~'**-^^+i~^^''^^'^ = a"^"'/"^. Choosing 
7 > 2/k, we can guarantee that the term c\loga\y/JIp{zk) dominates. 

Since the noise acting on the system is non-degenerate, we may assume the existence 
of constants D. ^ D_ > such that 



< Fr{f, (fi, z)Fr{f, (fi, zf ^ D+ . 



(D.52) 
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Thus the variance of M^^^_^ is bounded above by 



Zk + l 



2KZk 



2fe + l 



2KZk 



(D.53) 



A Bernstein-type estimate (cf. Lemma D.8 in Appendix D.4) provides the bound 

P{M^^^^^ <-x] ^ e-^'/2^+ . (D.54) 



Using (D.53) and (D.54) in (D.51) shows that we may assume Pk < 1/2, and the result 
follows from (D.47) and the definition of A^. Note that = log 2/7 < (log 2/2) /t. □ 



D.3 Laplace Transforms 

In order to combine the results from the two previous subsections, we will use a lemma 
based on Laplace transforms. In the following we let {xt}t^o be a time-homogeneous 
M''-valued Markov process with continuous sample paths. All subsets A C M*^ considered 
below are assumed to have smooth boundary, and to be such that the first-exit time 
ta = inf{t ^ 0: Xt ^ A} is almost surely finite. The Laplace transform of ta is the 
non-decreasing function 



[e^^^] = 1 + A / P^{rA > t} e^* dt G [0, 00] 
Jo 



(D.55) 



Note that W[e^^^] ^ 1 for all A ^ 0. Thus there exists a Aq ^ such that W[e^^^] < 00 
for all A < Aq. 

Lemma D.5. Choose nested hounded open sets 5i C ^2 C D C M'^. Let 7Z = V\ Si and 

consider the Laplace transforms 

GviX) = supE^[e^^^] , 

xeSi 

Gs{X) = supE^[e^^^2] , 
G^(A) = sup E^'[e^"'^] , 
g(A) = sup E^'[l|,^.<,^|e^^'^] . (D.56) 

Let A be such that Gs{\) and Gti{\) are finite, and assume that Q{X)Gs{\) < 1- Then 
G'x)(A) is also finite and satisfies 

Proof. For Xq E Si, one necessarily has ^ r-p, and thus the strong Markov property 
implies 



^xo ^gAri,j ^ ^xo gAr52 ^^r^^ ^gAri,j 
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Similarly, for x G 882-, since t-ji = t^c A tj) we have 



4f"Sf <f"'D} 



{-r-D<T5c} 



^g(A)G'2,(A) + G'7^(A) 



Now, ( |IX58| ) and <pl9\ i imply 

^^^(A) ^ ^^(A) [Q(A)G'2,(A) + Gn{\)] 

which yields the result. 



(D.59) 

(D.60) 

□ 



We will apply this lemma to sets Si = S{hi), S2 = S{h2) and V = 'D{k), where 
hi = Cj(j|log(T|, i = 1,2, with < Ci < C2. We introduce a new time variable t = z^, and 
let Xt be the time-homogeneous Markov process 

Proposition D.l yields a control of (^^(A) in the following way. The bound (D.6) 
translates in terms of the new process Xt as 

where Ci = C(z/)(/i2/cr)^'^ and Ai = K{u)/fi. It follows that 



(D.61) 



CiA 



POO 

Gs{X) ^ 1 + A / e^* sup P^°{r5, ^ t} dt ^ 1 + — 

holds for all A < Ai = k{u)/ ^. In a similar way. Proposition D.4 yields 

C2A 



A 



A, -A 



(D.62) 



(D.63) 



for all A < A2 = K,2/ ^A^oga\, where C2 = 2. It remains to estimate (5(A). Let us first show 
that Q{\) can be bounded in terms of (5(0). 



Lemma D.6. For all A < A2, one has 



1 - A/A2 



Proof. First note that for all T ^ 0, 



^ e^^ +A 



e^* dt . 



Plugging this into the definition of Q{X) yields 



(5(A) ^ sup 

x&dS2 



/oo 
^""{rn > t}e^* dt 



(D.64) 



(D.65) 



(D.66) 



The first term on the right-hand side is bounded by (5(0) e^^. The second one can be 

C2A 



estimated by Proposition |D.4[ yielding 

g(A) ^Q(0)e^^ + 



-(A2-A)r 



A. -A 



(D.67) 



Optimizing over T, we find that the optimal bound is obtained when e'*'^^ = (72/(5(0), 
which yields (D.64). □ 



50 



Finally, Q{0) can be estimated in a similar way as in the proof of Proposition D.4 
Proposition D.7. We have 



g(0) = sup P'jr^c <rv} 



K(c2-ci)^|log(7|2o/3(2;o)^/-D+ 



(D.68) 

Proof. Inequality (D.44) and the fact that e'^^^/^'^ / p{z) is increasing for sufficiently small z 



imply 



a 



{h2 - h)p{zo) + — 

71 



(D.69) 



for z ^ T-ji, where M° is the martingale introduced in (D.50). The variance of M° being 



bounded by D+p/2kzo (cf. (D.53)), the Gaussian tail estimate of Lemma D.8 allows to 
write 



P{rcc <zATn}=r\ inf (f, - hip{z)) ^ o] 

(C2 - C,f 



^ exp< —K 



\\oga\ 



K.{c2-ciY\\oga\zop{zo)'^ / Dj, 



(D.70) 



Note that the right-hand side of (D.70) does not depend on z. The result thus follows 
from taking the limit 2 — > 00. □ 



Proof of Theorem 6.4 Since A2 < Ai for sufficiently small cr, we set A = (1 — 9)X2 for a 
fixed < 6' < 1. For this A, we have 



C2 
6 



Ai — A2 



|loga| 



Furthermore Lemma D.6 and Proposition D.7 yield 







(D.71) 



(D.72) 



for some constant c > 0. Thus Lemma D.5 can be applied to show that G'd(A) is finite. 
Finally, by Markov's inequality. 



P^jr^, ^t} = P^'{e^"^ ^ e^*} ^ e-^*G'2,(A) 
which gives the theorem when translated back to the process (z- 



(D.73) 

□ 



D.4 A Gaussian Tail Estimate for Martingales 

Let Wz be an n-dimensional standard Brownian motion, and consider the martingale 



i=l 



(D.74) 



where g = {gi, . . . ,gn) takes values in M" and the process is assumed to be adapted 
to the filtration generated by Wz- We will assume that the integrand satisfies 



G^izf ^ giXz,z)giXz,zf ^ G+izf 



(D.75) 
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almost surely, for deterministic functions G±{z), and that the integrals 

V±(z)^ f G±(s)^ ds (D.76) 
Jo 

are finite. Note that under these conditions, is indeed a continuous martingale. 
Lemma D.8. For any x > 0, 

p| sup Ms>x \ ^ ^-xy2V+{z) _ 

Proof. Let 

[Ml^ [ g{Xs,s)g{X,,sf ds (D.78) 
Jo 

be the increasing process associated with M^. Then, for any 7 e M, the Doleans expo- 
nential 

g7M^-72[M]^/2 (D.79) 

is a martingale. It follows that 

pj sup > xl = pj sup e^^^ > e^4 

^ P| sup e^^^-^'[*^l^/2 ^ g7-7^[M]./2| 

^ p( sup > e^--^'^+(^)/4 



g7M,-72[M],/2 



(D.80) 



by Doob's submartingale inequality. Now the expectation in the last line is equal to 1, 
and the result follows by optimizing over 7, that is, choosing 7 = x/V+{z). □ 
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